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ABSTRACT 

The present invention relates to a computer-aided method to estimate the coil 
volumes using only two orthogonal angiographic projections and Noninvasive 
Vascular Elastography: Theoretical Framework. The invention relates to the 
validation on tissue mimicking phantoms of non-invasive elasticity imaging in small 
vessels. The Lagrarigian speckle tissue model estimator was used to assess the 
2-D-strain tensor, and the composite Von Mises elastograms were computed. The 
invention also relates to a model-based approach devoted to endovascular 
elastography (EVE). A modelrbased approach is proposed and the Lagrangian 
speckle model estimator is provided and the theoretical framework is validated. 
Moreover, the invention relates to the in vitro experiments using excised human 
carotid artery to validate the potential of endovascular elastography to detect hard 
plaques. 
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IMAGING APPARATUS AND METHODS 

FIELD OF THE INVENTION 

This invention relates generally to the field of ultrasound electrograpfry and more 
particularly, but not exclusively, to a rion-invasive ultrasound elastography system 
and a method for characterising the mechanical properties of superficial arteries. 

BACKGROUND OF THE INVENTION 

10 

In the early nineties, Ophir ef ai introduced elastography, which is defined as 
biological tissue elasticity imaging. Primary objectives of elastography were to : 
complement B-mode ultrasound as a screening, method to detect hard areas in the j 
breast Basically, the tissue under inspection is externally compressed and the ; 
displacement between pairs of pre- and post-compression radio frequency (RF) \ 
lines is estimated using cross-correlation analysis. The strain profile in the tissue is ; 
then determined from the gradient of the axial displacement field. 

1.1 Non-invasive vascular elastography (NIVE) 

Elastography has also found application in vessel wall characterization. However, ; 
20 vascular elastography is invasive; it is known in the literature as intravascular 
elastography or, sometimes, as endoyascular elastography (EVE). In EVE, the ; 
vascular tissue is compressed by applying a force from within the lumen. . Indeed, ! 
the compression can be induced by the. normal cardiac pulsation or by using a 
compliant intravascular angioplasty balloon. 

To validate the feasibility of EVE, phantom studies have been conducted. In most \ 
cases, tissue-mimicking phantoms with typical morphology and hardness topology 
synthesizing atherosclerotic vessels were constructed, de Korte et ai concluded 
that EVE may allow identifying hard and soft plaques independently of the 
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echogeneity contrast between the plaque and the vessel wall. The. potential of 
such an approach was emphasized by the fact that it provides information that 
may be unavailable from intravascular B-mode ultrasound (IVUS) alone. 

In vitro studies with excised vessels were also conducted for further validation of 
EVE feasibility. For instance, de Korte ef at computed elastograms from diseased 
human femoral and coronary arteries. They found different strain values between 
fibrous, fibro-fatty and fatty plaques, indicating the potential of EVE to distinguish 
different plaque constituents. Such results were compared with IVUS echograms 
10 and were corroborated with histology. They found that the elastograms were 
capable of demarcating regions within the « plaque representing differences in 
strain, whereas in IVUS echograms, these regions could not be discriminated. 
Using excised postmortem carotid arteries, similar results were recently observed. 

1 .2 Non-invasive vascular elastography (NIVE) 

So far, vascular elastography is invasive. Its clinical application is restricted to a 
complementary tool to assist IVUS echograms in preoperative lesion 
assessments and to plan endovascular. therapy. Our group recently, investigated 
the potential .of elastography to rion^nvasiveiy characterize mechanical properties 
of superficial arteries; the approach : was labeled as "non-invasive vascular 
20 elastography" (NIVE). 

Investigating the forward problem in NIVE, it was shown that motion parameters 
might be difficult to interpret; that is because tissue motion occurs radially within 
the vessel wall while the ultrasound beam, propagates axially. As a consequence 
of that, the elastograms are subjected to hardening and softening artifacts, which 
are to be counteracted. With NIVE, the Von Mises coefficient was propoised as a 
new parameter to circumvent such mechanical artifacts and to appropriately 
characterize the vessel wall The Lagrangian estimator was used, that is because 
it provides the full 2D-strain tensor necessary to compute the Von Mises 
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coefficient. The theoretical mode! was validated with.biomechanical simulations of \ 
the vascular wad properties. 

Simultaneously to the N1VE approach development, another group addressed the 
feasibility of strain imaging of internal deformation. A tissue-like gelatin elasticity- 
flow phantom was used to develop ultrasonic strain imaging for the detection of 
internal pulsatile deformations. Such an imaging technique also wqs applied in 
vivo to monitor deformation in tissues surrounding the normal brachial artery. 
Results suggested that the method may be feasible to detect regional artery 
10 elasticity changes caused by the formation of plaques and calcification. Although 
this last study investigated the -deformation of soft tissues surrounding an artery, 
the current study represents the first in vitro experimental results on the non- 
invasive characterization of smalt vesselmechanical properties. 

SUMMARY OF THE INVENTION : 

According to the present invention there is provided a method for characterising 
the mechanical properties of tissues in small vessejs using the Non-Invasive 
Ultrasound Elastography (NIVE), the method comprising the steps of : 
20 - Acquiring pre- and post- deformation images 

- Using the Lagrangten speckle model estimator to assess the 2-D stain 
tensor 

- Computing the Von Mi$es elastrograms 

- Mpdelising the tissue mechanical characteristics 

According to the . present invention there is also provided an apparatus for 
investigating the feasability of NI\>E for characterising tissues in small vessels. 

The feasibility of NIVE is investigated for the specific purpose of studying small 
30 vessels in humans and transgenic rodents with vascular diseases. Experiments 
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were conducted on tissue-mimicking phantoms. The equivalent vessel wall was \ 
made of polyvinyl alcohol cryogel (PVAO) and it was inserted into: a close-loop : 
pressurized device. 

The apparatus and the acquisition system are detailed in the section below. 

DETAILLED DESCRIPTION OF A PREFERRED EMBODIMENT OF THE 
INVENTION 
2.1 Experimental set-up 
10 Figure 1 gives an overview of the experimental set-up used to produce mechanical 
deformation of the PVA-G vessel-mimicking phantom, and the procedure utilized 
to collect RF. ultrasound raw data, allowing to compute vascular elastograms. A 
mixture of water-glycerol was circulated in the flow phantom. The height difference 
between the top (1) and bottom reservoirs (9) allowed to adjust the gravity driven 
flow rate and static pressure within the lumen of the phantom (3).. A peristaltic 
pump (10) was used to circulate the fluid from the bottom to the top reservoir. The 
flow rate (7) was measured with an electromagnetic flowmeter (Cliniflow II, model 
FM 701D), and the pressure (&) was monitored by a MDE Escort instrument 
(Uniflow pressure transducer from Baxter). As illustrated in Figure 1, the flow 
20 phantom was not directly connected to the tubing of the top reservoir to facilitate 
the small incremental pressure step adjustments necessary to obtain correlated 
deformation of the RF.signals within the PVA-C vessel wall. As described later on, 
the elastograms were computed from RF images obtained at two incremental 
pressure values. 

2.2 High frequency imaging system. ; (labels 2 and 4 on Figure 1) 

To allow • computing vascular elastograms of small vessels, a high resolution 
ultrasound , system with access i to the raw RF data was required. (Ultrasound 
biomicroscope, Visualsonics, ySAO, Toronto, Canada). This system was 
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equipped with a single element oscillating transducer with a central; frequency of • 
32 MHz, f-number = 2, diameter = 3 mm, focal length = 6 mm, and bandwidth at - ; 
6 dB of 110 .%. As recently reported, the axial and lateral resolutions of the \ 
instrument are 30 urn and 70 urn,' respectively. A transmit frequency of 40 MHz \ 
was used and the frame rate was 8 Hz. The bandpass filter was selected to : 
optimize the bandwidth.of the received echoes, and thus the axial resolution. 

2.3 Vascular phantom (label 3 on Figure 1) 

As shown in . Figure 2, the PVA-C vessel of the flow phantom was positioned 
between two watertight connectors, in a Plexiglas boxfilled With degassed water at . 
room temperature. Rubber brings were used to tight the PVA-C vessel onto 
Plexiglas tubes at both extremities. • 

As performed by others, the tissue-mimicking vessel was made of PVA-C 
(polyvinyl alcohol cryogel). This biogel solidifies and acquires its mechanical 
rigidity by increasing the number of freeze/thaw cycles. The freeze/thaw cycles 
modified the structure of the material by increasing the reticulation of fibers. It was 
shown that the elastic and acoustic properties of PVA-C are in the range of values 
found for. soft biological tissues. More specifically. Chu and Rutt demonstrated that 
the stress-strain relationship is very close to that of a pig aorta. 

The vessel-mimicking phantom investigated had a 1.5-mm lumen diameter, 1.5- 
mrn wall thickness, and 52-mm length. A 1.5% in weight of Sigmacell was added 
to the PVA-6 to provide acoustical: scatterers. The phantom simulated a double- 
layer vesselwall of approximately 0.75-mm thickness each. The stiffness of the 
inner wall was made softer by . controlling the number and temperature of the 
freezefthaw cycles. For this study, the numbers of cycles were set at 2 and 8 for 
the inner and the outer portions of the wall, respectively. 
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Figure 3 shows a schematic representation of both moulds that were used to ' 
produce the, doublerlayer vessel-mimicking, phantom. At a first instance, PVA-C j 
was poured between the first and second templates; that underwent 6 freeze/thaw « 
cycles to provide the external layer. At a second instance, liquid. PVA-C was ; 
poured between the second and third templates; that underwent 2 freeze/thaw j 
cycles to provide the double-layer vessel-mimicking phantom. 

: 

2.4 Acquisition protocol 

The ultrasound biomicroscope is a PC-based system (label 4 on Figure 1) ; 
providing, RF data that were pre-amplified (Panametrics, model 5900PR, label 5 on 
10 Figure 1) before. analog-to-digital conversion, After amplification, the signals were • 
digitized at 500 MHz with an 8-bit format by using a PC-based acquisition board . 
(Compuscope, Gage-Techmatron, model CS-8500, label 6 on Figure 1). 

2.5 Methods 

NIVE is based on the Lagrangian speckle model estimator (LSME) developed , 
already. The LSME is a 2D model-based estimator that allows computing the full 
2D-strain tensor. As illustrated in Figure 4, in NIVE, the observer's coordinate • 
system is the Cartesian (x,y)-plane while the motion coordinate system is the ; 
radial (r,<p)-plane. In such a situation, the parameters of an estimator are expected 
to be very difficult to interpret. The approach thus assumes a constant deformation . 
20 for small regions of interest (ROi), equivalent to W mn in Figure 4. If translation is « 
appropriately compensated for, this leads.to a simple linear transformation (LT) of i 
the ROI. The LT matrix relates the strain tensor s as follows: 
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The Von Mises (VM) coefficient was derived from the strain tensor s (Eq. 1) and it 
was proposed to counteract elastogram mechanical artifacts (such as hardening 
arid softening artifacts). The VM coefficient can be expressed as: 



\ - ^ 2 U +A 2 a -A H A 22 +|(A,2+ Aj,) 2 



(2) 



Resuits 

This section reports preliminary results for the double-layer phantom presented 
above. The pre-load static pressure was set at 3 mmHg (flow-rate of 70 ml/rnin), 
and pre-motion radio-frequency (RF) images were acquired. The static pressure 
was then increased at 7 mmHg, and post-motion RF images were recorded. 
Those RF images had a dimension of 8 mm x 8 mm, that is 5888 samples x 513 
RF lines. To apply the LSM'E, a measurement-window (W mn in Figure 4) of 272 jam 
x 312 urn (200 samples x 20 RF lines) was chosen. The motion was assessed for 
windows W™ with 85 % and 85 % axial and lateral overlaps, respectively. 

Pre-.and post-motion images were pre-processed to compensate for motion of the 
experimental apparatus, and then correlation coefficients were computed. It was 
then possible to examine the most relevant data among all experiments 
performed. Figures 5a and 5b show pre- and post-motion B- mode images; they 
are displayed for illustration purpose since computations were done on the RF 
images. Figure 5c presents a map of the correlation coefficient between the pre- 
and post-motion RF images; an average close to 0:6 clearly allows demarcating 
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the vessel wall from the lumen and the surrounding water environment Figure 5d 
shows the composite elastogram (4 in Eq. 2) as computed by the LSME, Because 
of the limited lateral resolution of the imaging system, motion estimates seem 
specifically reliable in the axial portion of §, delimited with dotted lines. For clarity, 
this axial part of the image is displayed in Figure 5e., showing clearly the inner 
layer (softer) with higher strain, values (close to 3 %), whereas. the* outer layer 
(harder) shows lower strain values (around 1%) as expected. This observation is 
more explicit in Figure 5f, where an average of 20 axial lines chosen in the middle 
of 5 is plotted. 

10 

This invention addressed the feasibility, of non-invasive vascular elastography for 
the purpose of investigating, small vessels in humans and transgenic rodents with 
vascular diseases. Experiments were performed in vitro on a double-layer vessel- 
mimicking phantom of 4.5-mm external diameter. The. Lagrangian speckle model 
estimator was used to assess th£ 2D-strairi tensor, and the composite Von Mises 
elastogram was computed. The two-layer vessel walls were clearly identifiable. 
The feasibility of NIVE for small vessel elasticity imaging was demonstrated. To 
our knowledge, these are the first experimental results reported on the non- 
invasive characterization of small vessel mechanical properties; arid this was 
20 made possible with the use of . a high resolution ultrasound system by 
Visualsonics. It may be noted, However, that the limited lateral resolution of the 
ultrasound system restricts the- strain contrast mainly in an axial region where 
tissue motion mostly runs parallel to the ultrasound beam. 
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FIELD OF THE INVENTION 

This invention relates generally to the field of endovasculaf ultrasound 
elastography (EVE) and more particularly, but not exclusively, to a model-based • 
approach to investigate EVE. ' 

BACKGROUND OF THE INVENTION 

Intravascular elastography, or equivalents endpvascular elastography (EVE), was • 
10 introduced in the late nineties as a new imaging modality that aims to outline, the ; 
elastiG properties of vessel walls. In EVE, the vascular tissue is compressed by ; 
applying a force from within the lumen. Indeed, the compression can be induced : 
by the normal cardiac pulsation or by using a compliant intravascular angioplasty : 
balloon. Phantom studies were conducted to validate the feasibility of EVE, In \ 
most cases, tissue-mimicking phantoms with typical morphology and hardness 
topology synthesizing atherosclerotic vessels were constructed, de! Korte et al. • 
(1997) demonstrated that EVE may allow identifying hard and soft plaques , 
independently of the echogeneity contrast between the plaque and the vessel wall. 
The potential of such an approach Was then emphasized by th§ fact that it : 
20 provides information that may toe unavailable from IVUS B-scans alone. At the 
same period, Soualmi et al. (1997) used finite element analysis to investigate : 
vessel wall elasticity images and. to demonstrate the feasibility of EVE. 

In vitro studies with excised vessels were also conducted to further validate EVE. * 

For Instance, de Korte et al, (1998, 2000a) computed elastograms from diseased i 
human femoral and coronary arteries. They found different strain values between 

fibrous, fibro-fatty and fatty plaques, indicating the potential of EVE to distinguish , 

different plaque constituents. Such results were compared with IVUS echograms ; 

and were corroborated with histology. They found that the elastograms were : 

30 capable of demarcating regions within the plaque representing differences in : 
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strain, whereas in IVUS B-scans v these regions could not be discriminated. Similar ; 
results were observed by Brusseau et ah (2001) using excised postmortem carotid \ 
arteries. In vitro EVE experimentation was also conducted by Wan et al. (2001). , 
They used an optical flow-based algorithm to .estimate the displacement field from \ 
B-mode data collected from porcine arteries; they computed the elasticity modulus • 
distributions by solving a mechanical inverse problem. While low spatial resolution 
of envelope data remains a . limitation, the method seemed encouraging because 
of the highest accessibility of B-mode compared to radio-frequency (RF) 
instrumentation. 

10 . ' 

However, in vivo, the position of the catheter in the lumen is generally off center ; 
and may move in response to the flow pulsatility; moreover, the lumen geometry is • 
often not perfectly circular. In such conditions, the ultrasound beam does not run ; 
parallel with tissue displacements, and appropriate coordinate systems are ■' 
required to model both the ultrasound propagation and the tissue motion. A 
consequence of this misalignment is that substantial decorrelation between the 
pre- and the posttissue-compression echoes is induced and must be 
compensated for. Ryan and Foster (1997) then proposed to use a 2D : : 
correlation-based speckle tracking method to compute vascular elastograms. 

20 Additionally, to circumvent catheter movement, Shapo et aL (1996a, 1996b) i 
suggested to compute the tissue, motion in the reference frame of the lumen's 
geometric center of the angioplasty, balloon. As this reference frame depends only 
on the balloon shape, it can remove artifacts associated with probe motion. Both of ; 
those phantom investigations tended to. demonstrate the feasibility of EVE in the ; 
presence of motion artifacts, and its potential to provide new diagnosis information ] 
that may help in the functional assessment of atherosclerosis. 

Another strategy was proposed by de Korte et al. (2000b) to minimize artifacts due { 

to catheter .motion; it consisted in. using pre- and post-motion images near ; 

3 0 end-diastole for a pressure differential of approximately 5 mmHg. The computed in \ 
vivo elastograms could detect an : area composed of hardened material, which was 
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corroborated with MJS B-scaris that revealed a large calcified ar^a* This last 
method may become a standard procedure . if one considers EVE without the use 
of an angioplasty balloon- 
So far, EVE appears a very attractive and promising tool to characterize the 
mechanical properties of vessel . walls; but, in return, it is potentially limited by 
signal decollation. This paper describes a new modef-rbased approach for strain 
computation in EVE. Assuming that speckle can be seen as a material property, 
tissue motion is investigated in the Lagrangian coordinate system (6r f equivalentiy, 
10 in the material coordinate system) instead of the Eulerian coordinate system. The 
Lagrangian speckle model estimator (LSME) (Maurice and Bertrand, 1999a) is 
then used to; assess tissue motion. The optical flow equations are derived and are 
integrated in; the LSME process to provide accurate tissue motion . estimates. The 
proposed method is implemented using the Levenberg-Marquardt algorithm. The 
theoretical model is validated with biomechanical simulations of the vascular wall 
properties. 



2 0 SUMMARY OF THE INVENTION 

According to the present invention there is provided a model-based method 
devoted to endovascuter Cartography, the method comprising the steps of : 

A) peif orming a modeling in performing the sub-steps of : 
A1) modeling the polar static-image formation; 

A2) modeling the tissue-motion; 

A3) modeling the polar dyhamic^mage-formation; 

B) estimating tissue motion in performing the sub-steps of : 
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B1 ) assuming speckle is a material property; 

B2) providing LSI (Lagrangiari speckle images); 

•* 

B3) using the Lagrangian Speckle Model Estimator (LSME) to assess 
tissue motion; 

B4) solving a non linear minimization problem; 
B5) using the opticalf low equations; 

thereby providing an accurate tissue motion estimate for providing an 
elastogram of the tissue/ 

According to the present invention, there is also provided an apparatus for 
1 0 characterisinjg the mechanical properties of tissues, the apparatus comprising : 

- a CViS ultrasound scanner; 

- a 30-40 MHz mechanical: rotating single element transducer; 

- an oscilloscope; 

- a pressuring system; arid 

- a computer. 

Intravascular ultrasound (IVUS) is known to be the reference tool for preoperative 
vessel lesion assessments and for endbvascular therapy planning. Nevertheless, 
IVUS echograms only provide subjective information about, vessel wall lesions. 
Since changes in the vascular tissue stiffness are characteristic of vessel 
20 pathologies, catheter-based endovascular ultrasound elastography (EVE) was 
proposed to qutline the. elastic properties of vessel walls. In EVE, the catheter is 
subjected to. cyclic movements in. the lumen because of the pulsatile blood flow 
motion; additionally, its position. may be off center, and the lumen geometry is 
generally not circular. As a result of those particularities, the tissue displacements 
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and . the ultrasound beam may misalign and cause substantial decorrefation : 
between the pre- and the posMissu^motion signals. Accordingly EVE does \ 
require robust motion estimators. In this paper, a model-based approach is 
proposed to investigate EVE; it considers the speckle as a material property. Such ; 
an assumption leads to the derivation of the optical flow equations, which are I 
suitably combined with the Lagrangian speckle model estimator to provide \ 
accurate tissue motion assessments. The theoretical framework is validated with 
simulated RF.data of vessel wall, pathologies. Results show the potential of the 
method to dissociate atherosclerotic plaques and healthy vascular tissue. 

The method and the apparatus will be detailed in the section below. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 6.- Schematic illustration of the image acquisition process in EVE. The ; 
transducer is placed at the tip of the catheter and cross-sectional imaging of the • 
vessel is generated by sequentially sweeping the ultrasound beam in a 360° '. 
angle. In this ideal situation, the ultrasound beam runs parallel with the vascular 
tissue motion, i.e. in the (r;<p) coordinate system. i 

20 

Figure 7.- Image-formatibn model for a 20 MHz polar scan system; a) shows the 
beam profile as a function of depth; b) presents the simulated polar B-mqde image 
for an homogenous vessel section; c) is the IVUS simulated image. 

Figure 8.- Schematic representation of an Ideal" plaque. The Young's, modulus for ; 

the normal vascular tissue, was 80 kPa, while the plaque, (three times stiffer) was I 

set at 240 kPa. To emulate boundary conditions as provided by the surrounding ; 
environment, the Young's modulus for the surrounding tissue was set at 1 000 kPa. 

- * ' 

30 Figure 9,r a) In vivo IVUS cross-sectional image of a coronary plaque; b) 2D finite, j 
element mesh of the unloaded real geometry with spatial distribution pf the plaque ; 
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constituents. The Young's, modulus for the healthy vascular tissue (or.adventicia & 
media) was 80 kPa, while the dense fibrosis (three times stiffer) was set at 240 
kPa, and the cellular fibrosis at 24 KPa (ten times softer than the dense fibrosis). 

Figure 10.- Schematic implementation of the polar dynamic imaige-formation 
model. 

Figure 11.- a) Theoretical . radial/strain elastogram for the idealized plaque; b) 
theoretical radial strain distributions along the vertical and horizontal lines 
10 specified on (a); c) radial strain elastogram as computed with the LSME; d) LSME 
radial strain distributions along the vertical and horizontal lines specified on (c). 

Figure 12.~ a) Strain-decay-Gompensated LSME elastogram, showing substantial 
contrast improvement between hard and soft materials; b) vertical 1 D iplot from the 
elastogram showing a contrast ratio close to .3 between the plaque and the normal 
vascular tissue, as it can be expected; c) horizontal 1D plot from the elastogram, 
showing effective strain, decay compensation, and thus a substantial improvement 
of the contrast ratio. 

20 Figure 13.- a) Theoretical radial strain elastogram for thie real plaque, showing 
very complex strain patterns; b) and c) show vertical and horizontal ID plots from 
the elastogram, respectively. Strain decay is specifically observed at the inner 
portion of the* vessel wall. 

Figure 14.- a) Radial strain elastogram as computed with the LSME for the real 
plaque; b) and c) vertical • and horizontal 1D plots from the elastogram, 
respectively.; Because of strain decay, there is not st clear demarcation between 
cellular and dense fibroses, specifically in c). 

30 Figure 15.- a) Strain-cjecay-compeosated LSME elastogram for the real plaque, 
showing a substantial contrast .improvement; b) and c) vertical and horizontal 1D 
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plots from the elastogram showing morel effective contrast ratio between dense \ 
and cellular fibroses, after strain decay compensation. 

DETAILLED DESCRIPTION OF A PREFERRED EMBODIMENT OF THE 
INVENTION 

The forward problem in EVE is addressed; it is followed by the derivation of the [ 
tissue motion estimator, and by the biomechanical simulations of the vessel walls. ] 
Results are presented in the following section, followed by the discussion of the < 
10 results and the conclusions and perspectives to this work. 

METHODOLOGY 

Endovascular elastography (EVE) , is a catheter-based modality, which gives 
insights about mechanical properties of the vessel wall. Following the example of 
1VUS, the transducer is placed at the tip of the catheter and cross-sectional 
imaging of a vessel is generated by sequentially sweeping the ultrasound beam in ; 
a 360° angle. Mechanical parameters (radial strain, in this case) are estimated : 
from analyzing the kinematics of the vascular tissue during the normal cardiac 
cycle (or, in some instances* in response to an angioplasty-balloon push). Figure 6 < 
2 o presents a schematic illustration of such a. process. 

The forward problem in EVE 

The polar static-image-formation model 

i 

The image-formation modei,. presented here, is for a rotating beam (single element ; 
or an array transducer system) and. is thus expressed in polar coordinates. It is 
based on previous works by . Bamber and Dickinson (1.980). Such . a - model was \ 
used by Meunier and .Bertrand (1995) to study speckle dynamics; ;it was also ; 
considered to investigate Speckle motion artifacts (KaHel and Bertrand, 1994; 
Maurice and Bertrand, 1999b). Under assumptions such as space-invariance of 
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the imaging \ system, and plane: strain, conditions for, the motion (that is no 
transverse motion is involved), the following simple 2D model was used (Maurice 
and Bertrandi 1999b): 

I'(x,y) m h?(x,y)®z r (x,y), . ; (1) 

where r(x,y) is the RF image, fY(x,y) is the point-spread function (PSF) of the 
ultrasound system, ® is the 2D convolution operator, and z'(x,y) is the acoustic 
impedance function, which can be modeled as a white Gaussian noise (random 
10 distribution of : uncorrected scatterers within the region of interest). As : it will be 
seen further in this section, the notation u ' tt is only for convenience and it does not 
refer to any mathematical operator. 

Because the speckle dimension varies with depth for polar scan systems such as 
the one described in Fig. 6, Eq. 1: is valid only for small regions of interest (ROI). 
Accordingly, the linear image-formation model is formulated using the 
superposition integral, given by: ■ 

l(r,<p) - J{h(r ) 9,r^q)'Hr^cp')r'dry(p^ (2) 

20 

where r and <p are the radial .(d&pjfh) and angular coordinates, respectively; l(r,$), 
h(r,<p), and z(r,<p) are the polar RF image, polar PSF, and the acoustic impedance 
function mapped in polar coordinates, respectively. Furthermore; for a. polar scan 
system, h(r,<jj) can be considered angular-position invariant; therefore, Eq. 2 
becomes: 

l(r,<p) = JJii(r,r\<p-(p').z(r',cp')r'dr'dq)'. (3) 
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It is convenient to model the PSF as a . cosine modulated by a 2D Gaussian 
envelope; that is a simple approximation of the far field PSF. The mathematical 
formulation can be expressed as: 



h(r,r',<|> 



2 °l 2 <*i 

v <p 



cos 27i-^f-(r~r') , (4) 



where o r is a pulse length parameter; cr<p p a^(r) is a beam-width parameter (the 
dependence in V is because of the radial-position variance of the PSF); f tr and c 
are the transducer frequency and the sound velocity in soft tissue, respectively. 
10 For simplicity, it is assumed that %q, (r) is a linear function of r; it is expressed as: 

Gy(r) = JL-a 0 , r^R L (5) 

with R L being the lumen radius and or 0 the beam width at r = R L , In such a 
situation, the beam forms a divergent cone with a section that linearly increases 
with depth through the vessel wait it is important to notice that Eqs 4 and 5 define 
a very simple approximation of the PSF in the far-field. In practice, the near field 
beam profile is more complex. However, it can be assumed that the transducer is 
positioned near the middle of the lumen, so that the vessel wall is not in the near 
20 field 1) . Additionally to the divergence of the ultrasound beam profile, the curvature 
of the wave-fronts also raises the complexity of the polar scan model; this will be 
argued in .the discussion. 

In theory, cr v (i) is a continuous function of r. However, for simplification it will be 
considered as a piecewise constant function,, the n partitioning regions being 



,} Whereas the lumen is in practice in the near, field, it is here simulated as being in the far field (with a beam 
width a^r) - cr 0 ). this simplifying hypothesis does not alter the generality of the model, since tissue motion 
. assessment is not; relevant for the lumen. ' t * 
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defined by the boundaries 0 = Co < rj< T2< ... < r n . For convenience, each partition 
will be denoted as Pk, that is: 

Pk=] r k-i> r k3> k=0,... a n (6) 
A discrete approximation of cr 9 (r) is then defined as foilows: 



10 where is the mean radial distance of the interval P k . Now, assuming that the 
PSF is locally depth-invariant, i.e; on each P k , the linear model of Eq. 3 can be 
approximated by the following convolution form: 



I(r 3 <p) * ElCPk.^^SWPfc^JiSzCP^q))] (8) 
k k 



with: 



20 While the continuity of l<r f cp) at the boundaries of the partitions can be questioned, 
Eq. 8 presents a. simple approximation, model of a polar scan system. Figure 7 
illustrates the; implementation of such a model for a 20 MHz transducer with a 6.0% 
bandwidth at -3 dB and a besm width (width at half maximum * 2.35xc D ) of 0.1 
mm. it is important to remember that z(r,cp) is the acoustic impedance function 
mapped in polar coordinates. z(rjcp) is assumed, to be a continuum, which can be 
modeled as a normal statistical process. For the purpose of simulations presented 
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below, z(r,<p) was obtained simply :by generating a 2D, normally distributed random 
field. The lumen and the tissue, surrounding the vessel were assumed to be 
respectively 2.5 and 1 .67 times le^sechoic than the wall. 



Figure 7a shows the beam profile as a function of depth. The beam width 
increases linearly, being minimal in the lumen and maximal in the surrounding 
tissue; the partition number h was set at 5. Figure 7b presents a simulated polar B- 
mode image for an homogenous vessel section, whereas Figure 7c gives the 
equivalent.lVUS image (in Cartesian coordinates). 

10 The tissue-motion model 

For a small ROl, tissue motion car* be approximated by an affine transformation; 
this can be expressed in Cartesian coordinates as: 



["p(x 9 y 9 t)l -j 
[q(x,y,t)J 




(9) 



20 



where Q f is a function of time t (9$)). Equation 9 is the result of a translation (vector . 
[T r ]) and of a linear geometrical transformation of coordinates (matrix [LT]); it can 
also be seen as trajectories that! describe a tissue motion in a region of constant 
strain (Maurice and Bertrand, 1999a). Strain is usually defined in terms of the 
gradient of a displacement field; since p(x,y,t) and q(x,y,t) represent the new 
position of a point (x,y), the (u x? u y ) components of the displacement vector in the 
(x,y) coordinate system are giveri by: 



M m [pO^O-xI re, "I 
Fe: 2 -i 93. T 



(10) 



with; A 
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In the above equation, A can be defined as the Cartesian deformation, matrix. The 
6(j, which are the components, of the strain tensor s, are expressed in terms of, the 
A(| components as: 



s e (t) = ^Ai(t)+Aj|(4 



(11) 



Furthermore, the radial and tangential components of the displacement vector (u r , 
u^) in the (r, <p) coordinate system are, respectively, given as: 



Ujp -arctan 



Uy 



, with u x su x (x(r,<p) 3 y(r,q>)) and u y su y (k(r,<p),y(r,q>)) 



(12) 



The polar deformation matrix, labeled as §(t), then can be derived from Eqs 12 and 
10; it is given, as: 



«t)- 



3u v : 


3u p " 






da 
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fir 
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3u r ' 






9u r . 
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fir 
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5tt x : 
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9u r :' 
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9u y 




_d<p 


dr 









"dx 


dx" 






dcp 
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~9y~_ 




dr 



with.: 



x^rcoscp. 
y«rsiriq> 



(13) 



The map of the ^ (= s r r) component is known in EVE as the radial elastogram; % n 
can be expressed as a function, of jthe Cartesian deformation matrix (A) as: 
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0r ~ 



9Uj 
9u, 



9u v 



coscp 
siacp 



(14) 



For a small pressure from within; the vessel lumen, £ rr is expected to provide a 
cartography of relative tissue stiffness inside the vessel wall. 

The polar dynamic-image-formatton model 

The 2D polar dynamic-image-formation model for an in-plane tissue motion is now 
derived. It is worth to remember that z(r,cp) is. a map of the acoustic impedance 
function z'(x.y) in polar coordinates. That can mathematically be expressed as: 



o 



z(r,cp)=z'(x,y) 



x=rcos(p 
y-rsm<p 



(15) 



The affin© transformation on z'(x.y) can be set by only changing; the (x,y) ; 
coordinates. Without lost of generality, it is assumed that translation is absent or is 
appropriately^ compensated for, and [T r ] in Eq. 9 can thus be neglected. The 
compensation for translation can be done using correlation techniques; such a 
processing is known, in the literature, as companding (Chaturvedi et al., 1998a 
and 1998b). It is also interesting to notjce.that impressing [LT] on the tissue to 
simulate motion requires to compute the inverse transformation [LT 1 ] on the 
coordinates. Hence, for an (r.cp) in-plane motion, the 2D RF polar dynamic-image- 
formation model at time t becomes:; 



I(r,q>»t) a tl(P k ,9,t)=|:[h(P^9)®^(?k»^ 



with z LT; ..(Pk><p)= z LT -'M -^LT'^y^r 



(16) 



lreP k 



: coscp 
= rsin(p 
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In Eq. 16, z'^, (x,y) indicates a change iry coordinates for the function z'(x,y); 
that change involves the 2 x 2 matrix [LT 1 ], Similarly, z LT -,(P k ,<p) indicates a 

change in coordinates for the function z(P k> <p) where the 2 x 2 matrixi [LTp 1 ] is 

involved, implicitly, this means that [LT] (as well ajs [LT P ]) is invisible. This 
assumption is valid for incompressible continuum . 

Tissue motion: estimation 
Lagrangian speckle image (LSI) ; - 

Vascular tissue is subjected to complex movement. For instance, it may rotate, 
shear, stretch or compress in .the measurement plane. As a result of * such 
10 kinematics, the relative positions of the speckle change; this sets a fundamental 
limitation to : correlation-based tissue motion estimators. However, it was ■ 
demonstrated: that the effects of such position changes can be compensated for in 
the process of motion assessment by pre-processing the signals using temporal 
stretching (Ophir et aL, 1999) br companding (Chaturyedi et al M 1998a and 
1998b). Interestingly, the Lagrangian -speckle image (LSI) was introduced to 
describe the ultrasound signals compensated for tissue motion (Maurice and 
Bertrand, 1999a). For example, in the dynamic image-formation model of Eq„ 16, 

tissue motion is induced by applying the linear transformation matrix [LT^ 1 ] to 
z(r,<p); hence, the motion-compensated image, said the LSI, is obtained through : 
20 applying [LT,-] (the inverse of [LT^ 1 ]) to a post-motion RF image at time t given by . 
l(r,<p,t). For the polar dynamic image-formation model given by Eq. 16, the LSI 
(noted as lu g ) is then expressed as: 
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£h(P k ^)®z L ^,(P k ,9) 



1t p 



(17) 



In the above equation, [LT p is the : determinant of the matrix [LT P ]. 



The Lagrangia;n speckle model and the minimization problem 

As described in Maurice and Bertrand (199i9a)^ the motion compensated RF image 
is labeled as the Lagrangian Speckle image, because the correction for tissue 
displacement directly involves the! Lagrangian description of motion. For instance, 
it is expressed in Eq. 17 that the LSI brings back material points to the positions 
where they originally stood. Accordingly, a convenient model to formulate the LSI 
10 can be given, as: 



1(1,9,0)= I Lag (r,(p,t)+9i(r,(p,t) - [l(r,cp,t)j LT + 9?(r,<j>,t) . 



(18) 



where 9*(r,cp,t> can be seen as an error term.. The mathematical model for a tissue 
motion estimiator then can be formulated as: 



MEN 
LT„ 



l(r,<|>,0)^[l(r,<p,t)3 LXp f m . ^Ijfc^-Iuifcftt)! 



(19) 



The minimum is obtained using the appropriate [LT P ]. It is worth to remember that 
20 [LT P ] is a linear transformation matrix; it maps the Cartesian trajectories (Eq. 9) in 
a polar coordinate system; However, for a small ROI (Ar, A<p) that is far from the 
lumen center, motion equivalent^ can be Investigated using either polar or 



CA 02457X71 2004-02-09 



24. 

coordinate system. In other words, the following approximation can be done to 
compute the eJastogram- 

£*LT-I (20) 

where I is the 2D-identity matrix. ; 

The Levenberg-Marquardt nonlinear minimization 

Several gradient-based methods exist to numerically solve minimization problems 
as given by Eq. 19. .For instance,; the Gauss-Newton proceeds iteratively towards 
the minimum , using, a linear approximation of 9?(r,(p,t). Because !Lag(r,<p f t) and 
5R(r f q> r t) are implicit functions of ft; (Eq. 9), and without lost of generality, let's 
rewrite i Lag (r,<p,t) as Iu g (§} and SR(r,q> t t) as 3*(e), respectively; 9 is the 
vectorization of. [LT P J. At the k m iteration, one will have: 

with I 0 =l(r,(p,0) 

In Eq. 21, Q k ~~ l \$ computed at ^iteration; k-1, Aef is the . increment of the i lh 
component of .6 , IuelG; is the.Lagrangian image at iteration k-1, and are 

6X7 38; 

partial derivative of tag with respect to each component of 9. A more compact 
formulation of Eq. 21 is given as: . 
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Lag 



with J k-i 5=5 



39t 



Lag 
39 1 



^Lag 
36 ^ 



81 



•Lag 



99,. 



Lag 
58 2 



ai 



Lag 



.99 r 



(22) 



h T T (S'-'l and * are vegans of .0. W and B.Wotively; 
where I 0 .k,ag\° . 

[lt J i5th e ^ matd* and *««°*° ' 

Lallan *U (3' i-oW M- '« - - *» ^ FW ^ . 
a model, the leaawquare error solution is given as: 



(23) 



10 v*er* tne subscript T designates *e transpose operator. When **I M * not 
The LUA regularized inversion was implemented as. 

(24) 



. onri 1 is the identity matrix, it is interesting to 

Know in« Gauss-Newton methods. Indeed, as -* o. 

20 features of the gradient and of the Gauss Nem 
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A6 is given by the Gauss^Newten method, that allows fast convergence in the 

vicinity of the solution; and as. A* r* «>, AG becomes parallel to the steepest- 
descent direction, that allows to converge even when the initial guess is outside 
the region of convergence for other methods. In th$ current study, the identity 
matrix was used as the initial guess: 

The optical flow equations and the; Jacobian matrix 

The optical flow equations,, or material derivatives, give a relationship between 
measures in Eulerian and Lagrangian coordinate systems, , respectively (Horn, 

1986). For instance, lu g ip k ~ l ) cart be seen as a function that describes a material 

10 property. Assuming that suqh a materia! property is preserved with motion, the 

total derivative of lugfp k ~ l } can be expressed as: 

d8i ' dr d9i 9q> dGy SO; 

(25) 

S9j dr 6B\ ! ' Ap ' dB s 1 




As discussed above, under the assumption of a small ROI (Ar, A<p) that is far from 
the lumen center, tissue motioh equivalent^ can be investigated using either the 
Cartesian or polar coordinates. Additionally, for small motion, the gradient of 
^ag^" 1 ) should not be . significantly different from the gradient ifc^o). These 
20 hypotheses, in conjunction with Eq. ;9, lead to: 
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Eq, 26 gives the full expression for the 6 components of the Jacobian matrix [J k ^] 
(Eq. 22). This is an interesting result, since using the optical flow equations to ■ 
compute the gradient appears more robust than the conventional finite difference : 
method. It is thus expected that Suph a combination of the optical flow equations 
with the L&M; minimization, algorithm may lead to more accurate tissue motion 
estimates. This last approach was; used in the current study. 

Biomechanical simulations of vessel wall kinematics 
Model design and image analysis ; 

The computational structural analysis was performed oh one idealized coronary j 
plaque (Figure 8), and on one typical composite plaque identified from an in vivo 
1VUS image. (Figure pa). For the idealized case, the lumen cross-sectional area . 
(LA) was 3,08 mm 2 , the external .elastic membrane cross-sectional area (EEMA) 
was 12,32 mm 2 , the plaque + medija crossTsectional area (P+MA=EEMA-LA) was . 
9,24 mm 2 , apd the plaque burden (PB%^100*(P+MA)/EEMA) was 75 % (Mrntz et : 
al., 2001). Equivalent^, for the real case,, the LA was 3,78 mm 2 , the EEMA was 
34,29 mm 2 , and the PB% was 89 % 

The major difficulty in computational structural analysis based on in vivo imaging is 
to determine the unloaded physiological configuration of the artery; i.e., the . 
configuration when the artery is. subjected to no external load. This: configuration 
has to be' known for finite. element £FE) simulations. To. obtain this unloaded state, ; 
adenosine triphosphate (ATP) . (Striadyne®, Wyeth France Laboratories); was : 
injected to the patient, as previously described by Ohayon et al (2001). All ; 
contours in the iVUS image were ^manually traced (Figure 9b); These contours are ! 
those of th0 lumen border- media, adventicia, and plaque components (dense . 
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fibrosis, mixed! fibrosis, and cellular fibrosis?). The adventicia contour was : added in 
the simulation ; and it had a mean thickness of 350 jjm (Rioufol et al., 1999), so as 
to take account of its protective roje: against any radial overstretching of the artery 
(Rachev, 1997). The various contours were digitized using the Un-Scan-lt® 
software (Silk Scientific, Inc., Orem, { UT). 

Material properties 1 

For the two models, the materials were considered . as quasi-incompressible 
(Poisson ratios 0=0.49) and isotropic with linear elastic properties. The Young's 
modulus for the healthy vascular tissue (or adventitia & media) was 80 kPa 
10 (Williamson et al., 2003), while the dense fibrosis (much stiffer) was set at 240 
kPa, and the defiular fibrosis (softer Jhan . the dense fibrosis) was chosen at 24 kPa 
(Ohayon et al M 2001; Treyve et : : al;, 2003); As to emulate boundary conditions 
provided by : the surrounding environment, these two vessel sections were 
assumed to be imbedded in a 1000 kPa Young's modulus tissue. 

Structural analysis 

FE computations were performed using the ANSYS 5.7® software (Ansys, Inc., 
Cannonsburg,; PA); Static simulations of. coronary plaque under loading blood 
pressure were performed oh the geometrical models previously described (Figs 8 
and 9b). Nod^tl displacements were* set to: zero on the external boundaries of the 
20 surrounding tissue. 

The various regions of the plaque components were then automatically meshed 
with triangular (6 nodes) and quadrangular (8 nodes) elements. The FE models 
were solved under the assumption of plane and of finite strains. The assumption of 
plane strain was made because aitfal stenosis dimensions were of at least the 
same order of magnitude as the : radial dimensions of the vessel. Moreover, the 
assumption of finite deformation was required as the strain maps showed values 
of 30% for physiological pressures (Loree et al., 1992; Cheng et al., 1993; Lee et 
al, 1993; Ohayon et al.; 2001; Williamson et al., 2003). The Newton-Raphson 
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iterative method with a residual rlod^l tolerance of 4x1 0" 4 N was used to. solve the 
FE models. The calculations wer^ performed with a number of elements close to \ 
7200. 

This computational structural FE analysis was used to perform the kinematics of ' 
the vascular tissue. The dynanhic imiage-formation model . (Eq. 16) was : 
implemented using the Matlab software. (The MathWorks Inc, MA, USA,:ver. 6.0). j 
The process, to simulate; polar, radio-frequency <RF) images is schematically 
presented in figure 10, for a homogeneous (pathology-free) vessel wall; it can be 
10 summarized as follows. It started &yi generating in Matlab a scattering function that ; 
simulated the acoustical characteristics of a transverse vascular , section in • 
Cartesian coordinates; that provided z ! (x»y). The axial and lateral displacement \ 
fields were computed with Ansys find were applied upon z^x.y) to perform motion 
and then to provide z' LT -t (x,y). . 

The next step consisted in mapping . z'(x,y) and z' Lr .i.(x*y) in a. polar coordinate ! 

system (r,<p) to provide z(r,<p) and;z Tr) {r,<p) (Eqs 15 and 16). Both polar-mapped . 

* p 

acoustic impedance, functions were .then convolved with the polar PSF to provide •; 
polar pre- and: post-tissue-mption RF images (Kr.^.O) and l(r f cp,t), given by Eq. 16). : 
20 Those images, were used as inputs to. the LSME (Eq. 19). 

In summary, the static and the. dynamic image-formation models associated with a ; 
polar scan were derived in section 1I.A, whereas the LSME was adapted for EVE 
in 1KB, In the Results section below, this new approach is validated using, the ; 
biomechanics! simulations of the; vessel wall kinematics (for the "id$a! w and the ; 
Realistic" plaque geometries) presented in Beetle . 

Results . 
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The idealized vessel measured about 3.8 mm in outer diameter, whereas the RF 
images extended to 4 mm x 4 mm; the re^il case vessel measured about 7 mm in 
outer diameter, whereas the BF: images extended to 8 mm x 8 mm. For the 
purpose of simulations, the. intraluminal pressure gradient for both vesselr-wall 
geometries described above was siet such as the dilation at the inner wall was 
around 7 %. The PSF characterized a 20 MHz central frequency transducer, as ■ 
described in the forward problem; (Section II.A). The LSME (described in section I 
H.B) was implemented to assess tissue motion; measurement-windows of 0.38 
mm x 0.40 mm and 0.77 mm x 0:80 mm, with 90 % axial and lateral overlaps, : 
10 were used for .the idealized and the- realistic cases, respectively. For more details 
concerning the definition of the measurement-window required with the LSME, the 
reader is referred to Fig. 1 of Maurice et at (2003). 

Investigation of the "ideal" piaque pathology 

Figure 11a presents the theoretical radial strain elastogram, computed for the ! 
"ideal" pathology case. The pteque can slightly be differentiated from the normal : 
vascular tissue, whereas a region of higher strain valyes is observed at the right 
portion of the inner vessel wali; Those "mechanical artifacts" are a direct ]• 
consequence of the well known strain decay phenomenon 2 * (Shapo et ah, 1996a). . 
For a more ^quantitative illustration, are presented in Figure 11b plots from the : 
20 theoretical elastogram for two orthogonal orientations along x and y. Indeed, the ; 
vertical plot ( — ) shows low contrast between the plaque and the normal vascular 
tissue, whereas the horizontal plot ( — ) clearly points, out the presence of strain 
decay. 

Figure 11c presents the radial strain elastogram as computed with the LSME, = 
using simulated RF images. As for the theoretical, elastogram in Figure .11a, the 
plaque is slightly distinguishable from the; normal vascular tissue and a region of f 
higher strain .values is also, observed. The plots of Figure 11d confirm . such : 



2) Radial strain z„ is proportional to l/rl This decreasing of e,,. with depth is usually defined as strain decay. 
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observations. Notice that lower ; strain values were computed in the LSME 
elastogram than in theory, specifically at the inner wall; that is due to the 
windowing process required to assess tissue motion with the LSME. 

■ • • * 

For the purpose of compensating for strain decay, the LSME radial sjtfain 
elastogram was post-processed.: Indeed, s rr was modulated with a function 
proportional to the square of the vessel radius 3> . In Figure 12 a js presented the , 
strain-decay-compensated LSME elastogram, showing substantial^ contrast i 
improvement. For instance, the axiql plot of Figure 1 2b shows an effective contrast : . 
ratio close to 3 between the plaque and the normal vascular tissue, as it can be ;■ 
expected; equivalent^, Figure 12c alsd : shows some valuable contrast . ratio 
improvement 

Investigation of a "realistic" vessel wail pathology 

Figure 13a presents the theoretical radial strain . elastogram, computed for the 
"realistic" pathology case, Interestingly, complex strain patterns are observed; 
nevertheless, different regions can.be identified. For instance, since the ratio of 
Young's moduli between the dense: and the cellular fibroses was arbitrarily set to 
10, both of those materials can. be; distinguished. Less contrast is seen between 
the cellular fibrosis and the . healthy vascular tissue because their Young's modulus 
contrast was set to 3. As illustrated^ with vertical, and horizontal 1D plots from the 
elastogram (Figs 13b and 13c, respectively), strong, strain decay is observed 
specifically at the inner portion of the vessel wall. 

Figure 14a presents the radial strain elastogram as computed with the LSME, 
using simulated RF images. As for the theoretical elastogram in Figure 13a, very 
complex strain patterns are observed. Moreover, the dense and the cellular 
fibrosis tissues can be identified: However, while. less prominent than in the "ideal" 

3) Since radial strain decreases proportionally to 1/r 2 , the compensation for strain decay consists in 
multiplying e„. by ia function proportional t6 r ? (Shapo et aL, 1996a). However, for non symmetric vessel 
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case study, strain decay, remains a significant factor to compensate for to improve 
image interpretation. This is illustrated in Figs 14b and 22c, where vertical and 
horizontal 1 D plots from the elastogram are presented. Whereas low strain values 
clearly indicate the presence of stiff -materials in Figure. 14b, this is not the case in 
Figure 14c. 

In Figure 15a is presented the strain-decay-compensated LSME elastogram, ; 
showing substantial contrast improvement. Now, the vertical plot {Figure 15b) as » 
well as the horizontal one (Figure 15c) show more effective contrast ratio between ; 
10 dense and cellular fibroses, and between cellular fibrosis and the normal. vascular .« 
tissue. Moreover, it is interesting to : notice the presence of moderate strain values ; 
(around 0.6 to Q.8 %) at the. extremities of the plots; this characterizes regions of 
healthy vascular tissue, namely the media and adventicia. 

Discussion 

Pathological conditions of vascular tissues often induce changes in the vessel wall ; 
elasticity. For instance, plaque deposit stiffens the vascular wall and then 
counteracts its dilation under systolic blood pressure. Hence, investigating .; 
mechanical and elastic properties of the arteries seem suitable, to appreciate the 
dynamics of the arterial wall and: its pathologies. In this paper, a model-based 
20 approach devoted to outline the elastic properties of the vessel wall with / 
endovascular elastography (EVE) was presented. Results obtained from numerical : 
simulations establish the potential of such a method to reliably assess very 
complex strain patterns. 

About the forward problem 

Regarding the forward problem (FlP) in EVE, a polar; static image-formation model < 
was introduced. Taking into account the spatial variation with depth associated to : 
polar scan systems, this jmage^formation model was formulated using the 

geometries, as it is the case In this study, this process may require to segment the vessel lumen as to provide : 
the vessel geometric center. 
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superposition integral. The radial variation was conveyed with the beam width, that j 
increases as a linear function, of depth; This is actually a relatively simplistic * 
approximation. In practice, the transducer point-spread function (PSF) is expected ; 
to be more complex. Indeed, the PSF geometry could be more curved than linear ! 
(Kaliel and Bertrand, 1994; Mauripfe and Bertrand, 1999b); such a curvature may . ; 
appear very important for tissue motion estimation, since it was demonstrated that : 
it would induce speckle motion artifacts (Kaliel and Bertrand, 1994; Maurice and 
Bertrand, 1999b), It would be worth, in the future, to address the polar dynamic : 
image-formation model with, a more complex ("realistic") PSF geometry, 

10 About the tissue motion estimation ! 

To assess tissue motion, the Lagrangian speckle model estimator (LSME) was ; 
used. The LSME is a 2D model-based estimator that allows computing the full 2D- 
strain tensor. In this paper, it was adapted for EVE investigations. While the full 2D : 
polar strain tensors were computed, only the radial strain parameters were { 
mapped to provide the elastograitis. This was motivated by the fact that tissue : 
motion, in EVE, is expected to run. parallel with the ultrasound beam. However, in \ 
real-life situations, this is not always the case. Inherent motions of the vessel wall ; 
that are due to the pulsatile blood flow, for instance, may misalign the ultrasound ; 
beam and the tissue motion.. While a 2D model partially compensates for such a \ 
20 decorrelation, the adapted Levenberg-Maqiiardt algorithm presented here also '. 
should help in circumventing this drawback. Indeed, the Jacobian matrix, as 
implemented with the optical flow equations (Eq, 26), provides a more robust 
gradient estimation method than the finite-difference method generally used. 

As discussed above, the. pre- and post-motion RF images are, in practice, ; 
subjected to motion artifacts, due to, the complexity of the beam geometry for polar \ 
scan systems. Accordingly an estimator, that allows counteracting such motion ; 
artifacts, could be required, For instance, the application of the Lagrangian filter, : 
as described in Maurice and Bertrand (1999a and 1997), could be a potential 
30 solution. Indeed, the Lagrangian. filter was proposed to take into account the PSF : 
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term in the motion estimation process; it thus potentially allows to compensate for ! 
artifacts introduced by instrumentation. We look forward for the implementation of 
the f iltered-version of the LSME for future clinical applications. 

Conclusion 

A new approach devoted to endovascular elastography (EVE) was presented. As i 
a first step, a polar static image r formation model was described; it took into ; 
account the spatial variation of the .PSF with depth associated to polar scan ; 
systems, The Lagrangian speckle model estimator (LSME) was then adapted for 
EVE investigations. Indeed, the LSME was formulated as a nonlinear minimization 
problem, for which an analytical formuiaticin of the Jacobian matrix was derived. [ 
The hypothesis behind this model-based approach was that speckle is a material \ 
property. While the full 2D. polar .strain tensor was assessed, only the. radial strain 
parameter was mapped to provide the elastogram. This was motivated by the fact ; 
that, in EVE, tissue motion is expected to run parallel with the ultrasound beam. To • 
validate the method, biomechanicaj simulations of the vessel wall kinematics were \ 
conducted. The results allow believing in the potential of this ; model-based 
approach to reliably assess very complex strain patterns; interestingly, it provides : 
mechanical property insights that are occlusive to IVUS B-scans. 
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FIELD OF THE INVENTION ■ • 

This invention relates generally . to the field of endovascular ultrasound 
elastography. (EVE) and tinore particularly, but not exclusively, to an apparatus and 
a method to detect hard plaques in the blood vessels, 

BACKGROUND OF THE INVENTION 

Atherosclerosis, disease of the intima layer of arteries, remains a major cause of 
. mortality . in western countries.. This pathology is characterized by a focal 

10 accumulation of lipids, complex carbohydrates, blood cells, fibrous tissues and 
calcified deposits, forming a plaque* that thickens and hardens the arterial wall A 
severe, complication of atherosclerosis is : thrombosis, consequently to plaque 
rupture or fissure, that might lead - t according to the event localization, to unstable 
angina, brain or myocardial infarction, and sudden ischemic death. Plaque rupture 
is a complicated mechanical process, correlated with plaque morphology, 
composition, mechanical properties and with the blood pressure and its long term 
repetitive cycle. Extracting information on the plaque, local mechanical properties 
and on the surrounding tissues may thus reveal relevant features about plaque 
vulnerability. Unfortunately no imaging modality, currently in clinical use, allows the 

20 access to these properties. 

So far, diagnosis and prognostic of atherosclerosis evolution mainly rest on plaque 
..morphology and vessel steno.se degree. These information can be accurately 
accessed with Intravascular Ultrasound (IVUS) imaging, since this modality 
provides high resolution cross-sectional images of arteries, Accurate quantitative 
analysis of the disease ..is thus easily performed by precise measurements of the 
lumen area, arterial dimensions and dimensions specific to the plaque. Moreover, 
IVUS permits a characterization of plaque components, but roughly, in terms of 
fatty, fibrous or calcified plaques and with possible misinterpretations. This makes 
30 IVUS, alone, insufficient to predict the plaque mechanical behavior. However, 
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elastic properties of vessel walls can be derived from radio-frequency (RF) 1VUS : 
images, by integrating elastographic processing methods. Indeed, endovascular ] 
ultrasound etestography (EVE) is an in-development imaging technique which : 
aims to outline elastic properties: Of vessel walls. Its. principle consists of acquiring . 
sequences of cross-sectional vessel ultrasound images, while the vasqular tissue \ 
is. compressed by applying a force from within the lumen; strain distribution is then 
estimated by tracking, within, the signals, the modifications induced by the stress ; 
application. In practice, in EVE, such a stress can be induced by the normal 
cardiac pulsation or by using a compliant intravascular angioplasty balioon. 
10 ; 
As primary investigations on EVE feasibility, phantom studies were conducted. In . 
general, tissue-mimicking phantoms . with typical morphology and hardness 
topology synthesizing atherosclerotic vessels were constructed. Namely, de Korte : 
et al. have demonstrated the potential of EVE to identify different vascular tissue i 
structures' (hard and soft plaques); independently of the echogeneity contrast with 
the healthy vessel wall. The potential of EVE was then emphasized by the fact that ; 
it provides information that may be unavailable. from IVUS alone. 

in vitro studies with excised vessels were also conducted to further validate EVE 
20 feasibility, de Korte et al. computed elastpgrams from diseased human femoral \ 
and. coronary arteries. They found different strain values between fibrous, fibro- : 
fatty and fatty plaques, indicating the potential of EVE to distinguish different > 
plaque constituents. Those results were compared with IVUS echograms and : : 
were corroborated with histology. The elastpgrams were found capable of •: 
demarcating regions within the plaque representing differences in strain, whereas \ 
in IVUS echograms, these regions could not be discriminated. Using excised l 
postmortem carotid arteries,, similar results were also observed by Brusseau et al. 

In vitro experimentation of EVE wlas also conducted by Wan et al. f using porcine ; 
30 arteries. In this case, an optical flow method algorithm was used to estimate the • 
displacement from B-mode data. The elasticity modulus distributions were \ 
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computed within the framework of the inverse-problem solution. White low spatial : 
resolution of envelope data remains a .limitation, the conclusion was the same: • 
EVE may allow identifying dffferent vascular tissue structures such as hard and : 
soft plaques. 

However, many difficulties arise from in vivo applications of EVE. For example, 
one of them concerns the fact that the, position of the catheter in the lumen is 
normally off center and not parallel to the vessel axis, and the geometry, of the : 
lumen is generally not circular. In such conditions, tissue displacements may be : 
10 misaligned with the ultrasound beam, introducing substantial decorrelation 
between the pre- and the posMissue-compression signals. Ryan and Foster ■ 
addressed this problem and proposed, to compensate for misalignment between 
tissue displacements and the ultrasound beam, the use of a two-dimensional : 
correlation-based speckle tracking method to compute vascular elastograms. 

Another potential drawback, associated with EVE in vivo applications, stems from 
the eventual cyclic catheter movement in the vessel lumen; that is due to the 
pulsatile blood flow motion. This may constitute source of signal decorrelation 
between pre* and post-tiss.ue-compression signals. To circumvent such a catheter • 

20 movement, Shapo et al. propqsedto compute the tissue motion in the reference 
frame of the lumen's geometric center of the angioplasty balloon. This was justified • 
by the fact that such a reference fratne depends only on the balloon shape; it thus ! 
removes artifacts associated with probe motion within the balloon. Another 
strategy, proposed by de Korte et al. to minimize artifacts due to catheter motion, :■ 
was to acquire, pre- and post-motion images near, end-diastole; that is for a : 
pressure differential of approximately 5 mmHg. Their approach was validated in [ 
vivo on patients referred to percutaneous coronary intervention.. However, only \ 
one elastbgram was computed; their results were corroborated with IVUS ; 
echograms, showing an area of hard material because of low strain values were : 

30 observed. 
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Furthermore, vascular tissue is heterogeneous, and its kinematics is complex. For 
instance, it may be subjected to rotation, shear and strain. As a consequence to 
that, changes in the relative position of the speckles within the measurement- : 
window and the changes in speckle morphology.induce a noise component, which \ 
increases the variance of the motion estimates. In addition, because of the force 
from within the lumen (from ..blood pressure excitation or from the compliant ; 
intravascular: angioplasty balloon), the vessel wall is also under the distension ■ 
stress; out-of-plane motion will then occur in the kinematics between the pre- and 
the post-motion images. This is another potential source of signal decorrelation ; 
1 0 limiting the reliability of tissue motion estimators in EVE. 

So far, EVE seems potentially limited by signal decorrelation. Recently, a model- • 
based approach with adaptation for EVE applications was proposed. It assumed \ 
that speckle is a material property; such an assumption allows considering the ; 
speckle dynamics in the Lagrangian coordinate system instead of the conventional : 
Eulerian coordinate system. From : that; the optical flow equations were derived, \ 
and were suitably combined with the Lagrangian speckle model estimator (LSME) 
to assess tissue motion. The approach, detailed in, combines the potential : 
features of the LSME with those of the optical flow equations. It is worth to 
20 remember that the LSME is a 2tp-tissue~motion estimator, that is implemented *• 
through an! iterative procedure using the regularized nonlinear minimization; 
method knqwn as the Leveriberg-Marquardt (L&M) algorithm. The L&M algorithm . 
is known for its robustness,, since it combines the best features of the gradient and 
of the Gauss-Newton methods. Indeed, the Gauss-Newton solution allows fast 
convergence in the vicinity of the solution, while the gradient method allows to ; 
converge even when the initial guess is outside of the region of convergence for . 
other methods, the optical flow equations are used to compute the Jacbbbian ; 
matrix that is required to implement the L&M algorithm. 



30 



CA. 02457171 2004-02-09 



39 i 

summary of the Invention ■ 

According to the present invention there is provided a method for detecting the ; 
hard plaque in the blood vessels, the method comprising the steps of : 

-collecting a sequence of RF images while incremental adjusting the : 
intraluminal static pressure; 

-characterising the the hard plaque from the normal vascular tissue. 

According*, to the present invention, there is also provided an apparatus for 
detecting the hard plaque in the blood vessels, the apparatus comprising : 

a CV1S ultrasound scanner; ; 

a 30 MHz mechanical; rotating single element transducer; 

an oscilloscope; 

a pressuring system; and 

a computer, : 

EndovasGular ultrasound elastography (EVE) was recently introduced to ; 
supplement; 1VUS echograms in: the assessment of vessel lesions and for : 
endovascular therapy planning. Indeed, changes in the vascular tissue, stiffness : 
are characteristic of vessel wall pathologies* and EVE appears as a very ; 
appropriate imaging technique to outline the elastic properties of vessel walls. • 
Recently, a model-based approach was proposed to assess tissue motion in EVE. 
It specifically consisted of a nonlinear minimization algorithm that was adapted to • 
speckle motion estimation- Regarding the theoretical framework, such an 1 
approach considered the speckle as a material property; this assumption, then 
leaded, to the derivation of the optical flow equations, which were suitably 



! 

Cfl. 02457X71 2004-02-09 



40 

combined with the Lagrangian speckle model estimator to provide the full 2D polar 
strain tensor. In this study, the proposed algorithm was validated in vitro using a 
fresh excised human. carotid artery. The experimental set-up consisted of a CVIS 
ultrasound scanner, working with a 30 MHz. mechanical rotating single element 
transducer, a digital oscilloscope; Lecroy 9374L and a pressuring system. A 
sequence of RF images was collected while incrementally adjusting the 
intraluminal static pressure steps. The results showed the potential of EVE to 
characterize and to distinguish atherosclerotic plaques from the normal vascular 
tissue. Namely, the geometry as well as some mechanical characteristics of the 
10 detected plaque were in; good agreement with histology. The results suggested 
that there may exist a range pf intraluminal pressures for which plaque 
detectability is optimal. 

This study consists in an in vitro validation of the adapted LSME. Indeed, a fresh 
excised human carotid artery was investigated, the experimental set-up consisted 
of a CVIS ultrasound scanner, working with a 30 MHz mechanical rotating single 
element transducer, a digital oscilloscope Lecroy 9374L and a pressurizing sytem. 
A sequence of 11 RF images was collected while incrementally adjusting the 
intraluminal static pressure steps. The elastbgrams were computed using the 
2 0 adapted LSME. Comparisons, with histology, allow believing in the potential of EVE 
to characterize and to distinguish atherosclerotic plaques from the healthy 
vascular tissue. As it will be seen, the results also suggested that there may exist 
a range of intraluminal pressures for which plaque detectability is optimal. 

The method and the apparatus will be detailed in the section below. 
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DETAILLED DESCRIPTION OF ; A PREFERRED EMBODIMENT OF THE } 
INVENTION 

There are presented the experimental set-up description, the data acquisition ■ 
protocol, and a summary of the Lagrangian strain estimation technique. Results on 
a fresh excised human carotid artery are presented after, while the discussion of • 
the results and drawing of conclusions and perspectives to this work follow. 

10 Materials and methods 

Experiments; were performed with | a fresh excised human carotid artery. This • 
section gives an overview of the specific devise that allowed both pressure 
variations inside the arterial lumen and RF data acquisition; it also summarizes the 
baselines of the method. • 

Experimental set-up description 

The experimental set-up was mainly made of a CVIS ultrasound scanner, working 
with a 30-40: MHz mechanical rotating single element, a digital oscilloscope 
LECROY 9374L and a self-made pressuring system (Figure 16). Artery extremities 
were fixed to two rigid sheaths ; by watertight connectors. The intravascular 
20 catheter was introduced through thje proximal sheath into the lumen of the artery, 
and then through the distal sheath. The distal sheath was closed with a clamp to 
insure wateftightness of the' system. Because of the sheath rigidity and of the 
system watqrtightness, injecting fluid inside the system will result in an increase of 
the pressure inside the arterial lumen. A syringe was then connected to the 
proximal sheath and the inner pressure was increased or decreased by manually 
varying the fluid volume (precision : AV- o;oi ml) inside the lumen. The probe was 
fixed approximately , at the centre of the arterial lumen tank to two guiding 
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elements. This has been performed in order to limit probe motion and accordingly 
reduce geometrical artefacts. 

Data acquisition 

In vitro experiments with the fresh carotid artery werfc performed at room 
temperature. , At each static pressure, . a scan of 256 angles was. performed. : 
Sampling of. the data was phase-synchronised, with the top image synchroniser ■ 
and the RF signal synchronisation (external output of the CVIS ultrasound ' 
scanner). The top image synchroniser allows the user to select an angular position ] 
from which the. acquisition started; it thus, permits the acquisition of sets of images 
angularly aligned. RF data were digitised at a 500 MHz sampling frequency in 8 j 
bits format, stored on a PCMCIA hard disc in the LeCroy oscilloscope and { 
processed off line. 

i 

Methods 

The tissue motion estimator, that was used to compute the elastograms, is : 
described in details in. Assuming a small region of interest. (ROI) and small tissue [ 
motion, it can mathematically be formulated as the following nonlinear \ 
minimization problem; 

lt p 

Where (r,<p) defines the image, coordinate system, while T gives the time. I(r,<p,0) 
is the pre-tissue-motion RF image, and lujfjfrt) is the Lagrangian speckle : 

image (LSI) at time T. So far, , it is worth to remember that the LSI yvas defined as ; 
a post-tissue T motion RF imagte that was compensated for tissue motion, as to 
achieve the best match possible with ifoqM)), The appellation "Lagrangian" refers 
to the Lagrangian description of motion. The minimum of Eq. 1 is obtained using 
the appropriate [LT P ], a linear transformation matrix. 
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So far it is important to remember that, as demonstrated in, .for a small ROI (Ar, 
Aq>) that is far from the lumen center, motion can equivaiently be investigated 
using either polar or coordinate system. In other words, the following 
approximation can be done: 



AsLT-I 



(2) 



where I is the 2D-identity matrix; [LT] is a linear transformation matrix; and A can 
be defined as/the Cartesian deformation tensor. Furthermore, it is known that for a 
small ROI, tissue motion can be approximated by an affine transformation; this can 
be expressed in Cartesian coordinates as: 
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where G| is a function of time t (8i(t)) f Eq. 3 is the result of a translation (vector fT r }) 
and of a linear geometrical transformation of coordinates (matrix [LT]). Eq. 3 can 
also be seen; as trajectories that describe a tissue motion in a region of constant 
strain. Assuming that (Ux,u y ) represent the displacement field in the (x,y) 
coordinate system, [LT] relates the strain tensor (e) through the following 
relationship: . 
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The radial strain then becomes equivalent to S22 (r A 2 2.= ©e - 1). the map of e 2 2 S 
distribution provides the radia! elastpgram.. 



Eq. 1 was numerically solved with an iterative procedure, that uses the regularized 
nonlinear minimization method, known as the Levenberg-Marquardt (L&M) 
algorithm. The optical flow equations wer$ used to . compute the Jacobbian matrix 
([J]) that is required to implement the L^W! algorithm; [J] was derived in; at the k m 
iteration, it can be expressed as: ; 
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0i are motion parameters, as . given by Eqs 3 and 4: [l k ] is a K m x n w .matrix, where 
"m H is the number of pixels in the ROI and "n", the 6 elements . of the affine 
transformation. To complete, it was demonstrated in that, at the (k-1) th iteration, 
the Jacobian matrix can be implemented as; 
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Eq. 6 gives the full expression, for the 6 components of the Jacobian matrix [J k ] 
(Eq. 5), The following section presents results Obtained with the proposed method, 
using data acquired from the excised human carotid artery. 



Results 
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Results are reported here for experiments performed with the excised carotid 
artery. This artery was characterized by a thin atherosclerotic plaque. A set of 1 1 
RF images were acquired for consecutive increasing physiologic fluid pressure 
levels. Figure 17 presents a histological section along with an IVUS image, 
Histology (Figure 17a) shows a very small atherosclerotic plaque (located at 3 
o'clock), that is only restricted to. a confined angular sector. The coloration with 
saffron haematoxylin-eosin revealed that the plaque contained cholesterol crystals . 
and inflammatory cells. Notice that the IVUS image does not allow to clearly 
differentiate the plaque from the healthy vascular tissue. 

10 

To implement the method, a 526 jim.-x 781 ixm. (200 samples x 20 RF lines) 
measurement-window, with 86 % axial and 90 % lateral overlaps, was Used. For 
the purpose of compensating for strain decay, the elastograms : were post- 
processed. In other words, they were modulated with a function proportional to the 
square of the vessel radius. Furthermore the elastograms were low-pass filtered, 
using a 600 p,m x 500 jam (6 x 6 pixels) kernel Gaussian filter, Figure 18 presents 
an overview of 10 elastograms that were computed. For instance, Figure 18a is for 
the lowest intraluminal pressure, .while Figure . 18] is for the highest one; indeed, 
maximum strain values close to .6 % are' observed in Figure 18a, whereas the 
20 maximum is close to 3 % in Figure 18j. To summarize, elastograms in Figs 18a 
and 18j respectively are the least representative, whereas those from Figure 18c 
to Figure 18e present very good plaque detectability, accuracy in plaque 
dimensions, and significant contrast between plaque and surrounding tissue. In 
other words, it does exist an optimal, range of. intraluminal pressures for which 
tissue motion estimation appears optimal. 

Figure 19 presents comparisons between histology (Figure 19a) and elastogram 
(Figure 19b). This elastogram (the same as Figure 18c) was chosen simply 
because it appears in the range of optimally described above, indeed, looking at 
30 the plaque geometry, the dimensions are very close to those observed in 
histology. Furthermore, from a biomechanics! point of view, stiffness ratio between 
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plaque and healthy vascular tissue is close, to 3; this is quite reasonable, regarding 
ranges of Young's moduli for srtfoerosclerotic plaques [xx kPa, yy kPa] and for 
normal vascular tissue [aa kPa, bb kPa] in the literature. 
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Non-invasive Vascular Elastography: Theoretical Framework 

Changes in vessel wall elasticity may . be indicative of vessel : 
pathologies. It is known, for example, that the presence of plaque stiffens the : 
vascular wall, and that the heterogeneity of its composition may lead to plaque 
rupture and .thrombosis. Another domain of application where ultrasound 

10 elastography may be of interest Is the study of vascular wall elasticity to predict the 
risk of aneurysmal tissue rupture. In this paper, this technology is introduced as an 
approach, to non-invasively characterize superficial arteries. In such a case, a ! 
linear array ultrasound transducer is applied on the skin over the region of interest, ] 
and the arterial tissue is dilated, by the normal cardiac pulsation. The elastograms, ; 
the equivalent elasticity images, are computed from the assessment of the . 
vascular tissue motion. Investigating the forward problem, it is shown that motion . 
parameters wight be difficult to Interpret; that is because tissue motion occurs 
radially within the vessel wall while the ultrasound beam propagates axially. As a 
consequence of that, the eiastograms are subjected to hardening and softening ! 

2 0 artefacts, . wh)ch are to be counteracted. In this paper, the Von Mises coefficient is ' 
proposed as, a new parameter to circumvent such mechanical artefacts and to \ 
appropriately characterize the vessel wall.: Regarding the motion assessment, the 
Lagrangian estimator was used; that is because it provides the full 2D-strain 
tensor necessary to compute the Von Mises coefficient. The theoretical model was 
validated with biomechanical simulations of the vascular wall properties. The 
results allow; believing in the potential of the method to differentiate hard plaques 
and lipid pools from normal vascular tissue- Potential in vivo implementation of 
non-invasive vascular elastography. to characterize abdominal aneurysms and 
superficial arteries such as the femoral and the carotid is discussed. 
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Index termsi Ultrasound elastography, Vascular wall, Mechanical properties, Non- 
invasive scanning, Mathematical modeling,; Vascular pathologies. 

INTRODUCTION 

Pathological conditions often induce changes in biological tissue stiffness. ; 
That is, for example, the basic hypothesis supporting palpation as a screening : 
method to detect hard tumors in the breast, prostate and other organs. However, \ 
in many instances, such an approach is impracticable when deep organs are 
10 considered. This gave rise to elasticity imaging, which aims to outline the elastic ; 
properties of biological soft tissues using ultrasound. 

In the early nineties, Ophir ef &L introduced elastography, which is defined • 
as biological; tissue elasticity imaging. Primary objectives of elastography were to = 
complement B-mode ultrasound, as a screening method to detect hard areas in the 
breast, and : to investigate prostate cancers: Such approaches are sometimes ■ 
referred to as computed palpation. Basically, the tissue under inspection is . 
externally compressed and the displacement between pairs of pre- and post- 
compression radio frequency . (RF) lines is estimated using cross-correlation • 
2 o analysis. The strain profile, in the tissue is then determined from the gradient of the ; 
axial displacement field. ... . 

Elastography has also found application in vessel wall characterization, i 

Until now, vascular elastography is invasive and is known in the literature as i 
intravascular ultrasound elastography or, sometimes, as endovascular 
elastography (EVE). In EVE, the vascular tissue is compressed by applying a force 

from within the lumen. Indeed, the. compression can be induced by. the normal ; 

cardiac pulsation or by using a compliant intravascular angioplasty balloon. Such a : 

balloon may also be used to stabilize the EVE catheter position in the lumen; This ] 
30 approach can be seen as a complement to B-mode intravascular ultrasound 

(IVUS). : 
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REPORTED WORKS ON EVE (RELEVANT TO THE CURRENT STUDY) 

* 

Simulations and phantom studies 

One of the first investigations on endovascular elastography was presented ; 
by Soualmi et ah in 1997, indeed, they used finite element analysis to better 
understand the vessel wall elasticity images. To validate the feasibility of EVE, 
phantom studies have been proposed. In most cases, tissue-mimicking phantoms 
with typical morphology and hardness topology synthesizing atherosclerotic : 
10 vessels were constructed* de Korte ieff al concluded that EVE may allow identifying 
hard and soft plaques independently, of the . echogeneity contrast between the 
plaque and the vessel wall. The potential of such an approach was emphasized by . 
the fact that it provides information that may be unavailable from IVUS alone. 

In y/Vo, the position of the catheter in the lumen is normally off center, and i 

the geometry of the lumen is generally not circular. In such . cases, tissue ; 

displacements may be misaligned with the ultrasound beam, introducing : 

. substantial decollation between the pre- .and the post-tissue-compression ; 
signals. As to prevent such limitations, Ryan and Foster proposed the use of a 

20 two-dimensipnal correlation-based speckle tracking method to. compute vascular ; 

elastograms. While questioning the limitations of the lateral tracking performance, ; 

they concluded their phantom studies claiming that EVE may provide a new ! 
spectrum of information to aid in the assessment of atherosclerotic lesions. 

Additionally, in vivo, the catheter may cyclically move in the lumen because ? 
of the pulsatile blood flow motion. This constitutes another potential source of ; 
decollation for the pre* and the post-tissue-compression signal treatment. To : 
counteract that, Shapo et a/, proposed to compute the tissue motion in the • 
reference frame of the lumen's geometric center. of the angioplasty balloon. .Such a • 
30 reference frame depends, only on the balloon shape, it thus removes artefacts : 
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associated with probe motion within the balloon. Their phantom investigations also \ 
tended to demonstrate the feasibility of EVE and its potential to differentiate [ 
between hard and soft plaques. 

More recent results on endbvascular phantom studies were obtained by 
Brusseau et aL. Instead of the traditional correlation-based methods, they used an 
iterative approach to compute a. scaling factor between pre- and post-compression ; 
signal segments. From that, the. strain was estimated. This adaptive method is i 
expected to . be accurate in a wider range of strains than the commonly used ; 
10 gradient-based methods, and it may prove better for investigating highly 
heterogeneous tissues; 

In vitro studies with excised vessels 

One more step to validate the feasibility of EVE is in vitro experiments, de ! 
Korte et aL computed elastograms from diseased human femoral arid coronary ; 
arteries. They found different strain values between fibrous, fibro-fatty and fatty ; 
plaques, indicating the potential of EVE to distinguish different plaque constituents. 
Such resultf were compared with IVUS echograms and corroborated with , 
20 histology. One of their principal findings was that the elastograms were capable of 
demarcating regions within the plaque representing differences in strain, whereas 
in IVUS echograms, these regions could not be discriminated. Using excised 
postmortem carotid arteries, similar results were recently observed by Brusseau et 
aL 

In vitro endovascular elastography also was conducted by Wan ef aL They 
used an algojrithm-basecl optical flow method to estimate the displacement from B- . 
mode data collected from porcine arteries. They also computed the elasticity ; 
modulus distributions within the framework of the inverse-problem solution = 
30 (estimation of the Young's modulus). While low spatial resolution of envelope data ; 
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remains a limitation, the method seemed encouraging because of the highest » 
accessibility of B-rriode compared to :RF instrumentation. 

In vivo studies 

Currently, in vivo applications of EVE are scarce, de Korte et a/, attempted 
to compute elastograms in vivo from patttents referred to percutaneous coronary 
intervention. Indeed, as pre-interverition IVUS assessment of the lesions was 
performed, RF data were acquired to compute elastograms. To minimize artefacts 
10 due to catheter motion, prer and post-tissue-motion images were acquired near 
end-diastole; that is for a pressure differential of approximately 5 mmHg. One 
elastogram was presented; it identified an area as being composed of hard 
material since low strain values were found. Suqh a finding was corroborated with 
the IVUS echograms that revealed a large calcified area. That result allows 
believing in the potential of EVE for follows-up of patients with vascular diseases. 

NON-INVASIVE VASCULAR ELASJOGRAPHY (NIVE) 

So far, vascular elastpgraphy is invasive. Its clinical application is thus 
20 restricted to a complementary tool to assist IVUS echograms in pre-operative 
lesion assessments and to plan endoyascular. therapy. It has also found 
application in in vitro studies to characterize vascular tissues. Nevertheless, 
elastography is a very attractive and promising approach to characterize the 
mechanical properties of vascular walls; In this paper, the feasibility of non- 
invasive vascular elastography is investigated. 

In the. next sections, the Lagrangian estimator and its implementation for 
vessel wall characterization are summarized. The full 2D-stra«n tensor s is 
provided, but the so-called elastogram is given by the strain component syy, which 
30 corresponds to the strain in the direction of the ultrasound propagation. 
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Investigating; the forward problem, it is shown in section 2 that such an elastpgram 
is subjected to hardening arid softening artefacts; that is explained by the fact that 
motion occurs radially within the. vessel wall while the ultrasound beam propagates : - 
axially. To circumvent such incoherence, the Von Mises coefficient is proposed as 
a new parameter to nori-irivasively characterize the vessel wail. Results from 
biomechanics simulations of walltissue displacement validate the potential of the ; 
approach to :accurately differentiate hard plaques and lipid pools from the normal 
vascular tissue. Potential in vivo implementation of non-invasive vascular 
elastography to characterize abdominal aneurysms and superficial arteries, such 
10 as the femoral and the carotid, is discussed in section 4. 

METHODOLOGY 

In conventional elastography, an axial compression is applied so that the I 
tissue motion, occurs in the ultrasound beam axis. In such a case, 1D estimator 
can provide reliable tissue motion assessment In some circumstances, 2D : 
companding :may be required to compensate for non-axial motion. In non-invasive ! 
vascular elastography (NIVE), the challenge stems from the fact that- motion 
20 occurs radially within the vessel wall section while the ultrasound beam 
propagates axially (linear atrray transducers are considered in the present study). ; 
Regarding that, 1D estimators are npt appropriate. In this paper, the Lagrangian 
speckle model estimator (LSME) Is.thus proposed to non-invasively characterize 
vascular tissues. Indeed, the LSME is a 2D model-based estimator, that allows • 
computing the full 2D-strain tensor. The approach relies on a tissue-motion model 
and on a dynamic image-formation model that are summarized in the next ; 
sections. 



THE TISSUE-EMOTION MODEL 

30 
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Subjected to blood flow pressure, the apparent motion of an artery is dilation; 
however, beeause of the boundary conditions imposed by the surrounding tissues 
and organs, the vascular wail itself is compressed. Such a compression induces a 
radial strain that is maximum at the inner wall and rninimum at.the outer wall; this 
observation is referred to, in the literature, as the strain decay. However, as 
illustrated in Figure 20, motion occurs radially in NIVE and is parallel to the 
ultrasound beam only at angles <p of .O 0 and 180° (note that cross-sectional ifrtages 
are considered here). Moreover, the vascular tissue is heterogeneous; it is thus 
expected to deform non-uniformly. Accordingly, proceeding to motion estimation 
requires subdividing the region under study into, several 2D sub-regions of interest 
(ROI). This is. shown in Figure 21, where the ROIs are represented by the 
measurement-windows Wm n - 

For small ROIs (W m n)v tissue motion can be approximated by the zerororder 
and first-order terms of a Taylor-series expansion; this can be expressed as: 



r P (x ?y ,t)i re, in 

.Lq(x,y,t)J „ [9 4 J (_< 



8, 8 3 

e 5 e 6J 



(1) 



Tr LT 



where Q t is a function of time t (9i(t)). Equation 1 defines an affine transformation, 
i.e., it is the. result of a translation (vector [T r ]) and of a linear geometrical 
transformation of coordinates (matrix [LT]). Equation 1 can also be seen as 
trajectories that describe a tissue motion , in a region of constant strain. Strain is 
usually defined in terms of the gradient of a displacement field; hence, as p(x,y,t) 
and q(x,y,t) represent the new'position of a point (x,y)» the (u,v) components of the 
displacement vector in.the (x,y) .system, are giv<3n by: 
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(2) 



with:. A= 



The ey, which are the components of the 
terms of the Q parameters as: 



?ii(t) - |KW+Aji"(tj3 



strain tensor e, can then be defined in 



(3) 



For an ultrasound Rulse. propagating in the y-axis direction, Cxx.and eyy 
represent respectively the lateral and the axial strains; « %x is the shear strain. 
10 Because conventional elastography assumes a constant stress distribution in the 
ROI, the axial: strain is inversely proportional to the Young's modulus. Accordingly 
Syy, known in the literature as the elastogram, can be seen as .a map of the tissue 
stiffness. In Qther words, elastography may allow identifying different tissue 
structures such as hard plaques and lipid pools, or heterogeneity in the 
aneurysmal vessel wall. 

THE DYNAMIC IMAGE-FORMATION MODEL 



It is assumed that the image formation can be modeled as a linear space- 
2 o invariant operation on a scattering function. If it is further assumed that the motion 
occurs in plane strain conditions (that is noHransverse motion is involved), then the 
following simple 2D model can.be used: 



;i(x/y) - H(x,y)®Z(x,y) 
<4) 
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where !(x,y)is the RF image, H(x,y) is the ppintspread function (PSF) of the 
ultrasound system, ® is the 2D convolution operator, and Z(x,y) is the acoustic 
impedance function, which can be modeled as a. white Gaussian noise (random 
distribution of scatterers. within the : ROlj. 

Let's now derive the 2D dynamic. image-formation model for an in-plane 
tissue motion. For small ROIs such as in Fig. 20, motion impressed on the tissue ( 
Z(x,y) ) can be assumed an ^ffine transformation (Eq, 1). This means that such a 
motion can be set by . only changing the (x,y) coordinates. Without lost of 
10 generality, it is assumed that the translation [Tr] is absent (no rotation of the blood 
vessel during pulsation) or is appropriately compensated for, and pan thus be 
neglected. The compensation for translation can be done using correlation 
techniques; such a processing is known, in the literature, as companding. It is also 
interesting to notice that impress [LT]. on the tissue requires to compute the 
inverse transformation [LT 1 ] on the coordinates. Hence, for an (x,y) in-plane 
motion, the 2D RF dynamic image-formation model is given by: 

l(x,y,t) - H(x s y)®Z Lr : t (x;y) (5) 

20 where Z LT .,(xjy) indicates a change in coordinates for the function Z(x,y); that 

change involves the 2 x 2 matrix.JLT* 1 ]. Implicitly, this means that [LT] is invertible. 
This assumption is valid for incompressible continuum. Indeed, in such a case, the 
determinant of [LT] (also known as the Jacobbian) is unity, so [LT] is nonsingular 
and invertible. 

Biomephanidal simulation ottissue motion 

As a first step, the Matlab software (The MathWorks Inc, MA, USA, ven 6.0) 
was used to simulate the simple case of a . homogeneous vessel wall. The process 
30 is schematically presented in Fig. 21 and pan be summarized as follows. It started 
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. by generating in Mattel a scattering Junction that simulated the acoustical 
characteristics of a transverse vascular section; that provided Z(x f y)> The 
analyticalisolution for a pressurized thick-wail cylindrical blood vessel embedded in ; 
an elastic coaxial cylindrical medium was then derived (Appendix I). From that, the : 
axial and lateral displacement fields were computed. Those displacement fields 
were applied upon Z(x,y) to perform motion and then to provide z LT _,(x,y). Finally, ; 

Z(x,y) and z LT .,(x s y) were respectively convolved with H(x,y) (the PSF) to produce ; 

the pre- and the post-tissuerrmotion RF images (l(x,y,0) and I&y.t), t * 0). As it is • 
seen in the next section, those two images constitute the inputs to the Lagrangian : 
10 speckle model estimator. 

. A pathologic (heterogeneous) vessel wall was also investigated, that is the 
case where the lumen. is narrowed by the presence of a hard plaque at the intima. 
The process to simulate the pre- and the post-tissue-motign RF images was 
similar to the. homogeneous case study, except that the Ansys finite^element 
modeling software (Ansys Inc., Canonsburg, PA, ver. 6.0) was used to describe 
the static mechanical behavior and the kinematics of the pressurized 
heterogeneous vessel. 

20 THE LAGRANGIAN SPECKLE MODEL ESTIMATOR (LSME) 

The Lagrangian speckle, estimator was introduced and widely discussed, it is a 
model-based, estimator that relies on the previous tissue-motion and dynamic 
image-formation models It can . be formulated as a nonlinear minimization 
problem. The LSME iteratively computes the [LTJ that allows the best match within 
a sequence of images; it thus provides the strain tensor through Eqs 1, 2 and 3, 
For simplicity and without lost of generality, let's consider only two images of the 
sequence: l(x,y,t) at t = 0, that is t(x f y,0) ■ lo(x,y) which is labeled as the pre- 
motion image; and I(x,y,t) at a. given time (ti * 0) such that l(x,y,ti) = h(x,y), which 
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corresponds, to the post-motion image. The estimator is schematically represented 
by the block-diagram of Fig, 22 and it can be expressed as: 

MIN Hx>y)-i^{x>yf : (6) 

with \|/ - [Tr;Lt(:)] 4 . Ii a g(x,y) is labeled as the Lagrangian speckle image; it can be 
expressed as Ili{x,y)]tr. that is the post-motion RF image ii{x,y) compensated for 
tissue motion. In Fig. 22, this is implemented, through an adapted version of the 
Levenberg-Marquardt (L&M) rninimlzatibn algorithm. . 

10 

Such an estimator provides many advantages. The major one is related to the 
fact that it allows computing the full- ZD-strain tensor (Eq; 3). As A*x (- exx) and A w 
(= syy), the divergence parameters,, provide information about tissue stiffness, the 
shear parameters, A*y and Ay*, may provide useful insights on the heterogeneous 
nature of the tissue under investigation. However, one should be aware that the 
accuracy of the lateral parameters. (A** and Axy) is potentially limited by the lateral 
resolution of ultrasound systems (the lateral resolution, that depends on the 
ultrasound beam characteristics, is generally lower than the axial resolution 
determined by the transducer bandwidth and the ultrasound system electronic 
20 properties). 

THE FORWARD PROBLEM: TISSUE MOTION ANALYSES 

For a : continuum, motion can be described in a Lagrangian coordinate 
system or in an Eulerian coordinate system. In the literature, the Eulerian 
coordinate system is sometimes referred to as the observer's coordinate system, 
whereas the Lagrangian coordinate system is known as the material coordinate 
system. The material coordinates allow to express each portion of the continuum 



4 This is the.Matlab notation for augmented vector {;) and matrix vectorisation (:), Hence, \\s is a 6 x 1 vector 
built from the 2 k i Tr vector and the 4 x 1 vectorisation 1 of LT. ' 
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as a function of time and position. The difference between these two coordinate 
systems is illustrated in Fig/ 23/ where the (x,y) constitutes the observer's 
coordinates and the (r,q)) defines the material coordinates. 

.With most imaging systems, such as ultrasonography, the observer's and 
the material coordinate systems are generally the. same; hence t most tissue 
motion estimators use, by definition, the observer's coordinate system. However, 
the material ■■ coordinates were presented as a suitable way to describe speckle 
dynamics, lt;was demonstrated that such a coordinate system leads to appropriate 
10 signal processing, which allows to counteract speckle decorrelation effects in 
tissue motion assessment In non-invasive vascular elastography, as illustrated in 
Fig. 20, the observer's coordinate system is the Cartesian (x,y)-plane. This system 
is different from the motion coordinate system that is in the radial (r.(p)-plane. In 
such a situation, the parameters of an estimator are expected to be very difficult to 
interpret. The challenge of noninvasive vascular elastography greatly concerns 
the interpretation of the estimated motion parameters to characterize the vascular 
tissue. In the next sections, such motion parameters are investigated and a 
proposition is given for a new tissue characterization, parameter. 

20 Motion analysis for a homogenous tissue 

To initiate the forward problem, a pathology-free situation was considered, 
that is the simplest case of a circular, axis-symmetric and homogeneous vessel 
section. To take into account the: constraints induced by the environmental tissues 
and organs, ft is hypothesized that the. vessel section is embedded in an infinite 
medium of higher Young's modulus. An exact solution of a pressurized thiok-wall 
cylindrical blood vessel of inner and outer Jradii R\ and Ro respectively, embedded 
in an elastic coaxial cylindrical medium of radius R e , is derived for our study in 
linear elasticity. The details are presented in Appendix i. In this study, it is 
30 assumed that the plane strain condition for the vessel wall applies and also that 
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the two media are incompressible and isotropic; Referring to Appendix I, the 
displacement gradient components (Eq. 2) are given by: 



A(x,y)= 



y 2 -x 2 



-2KL 



xy . 



-2K. 



x 2 -y a / 



xy 



K. 



(7) 



with: 



(8) 



Where w arid <2) respectively describe the vessel wail and the external medium; Pb 
1 0 defines the blood pressure and E is the Young's modulus. 



20 



Using; Matlab software* Eqs 7 and 8 were processed to simulate the 
dynamics of a homogeneous vessel section subjected to an intraluminal pressure 
gradient of 73 mmHg; that is close, to 6 % : intraluminal dilation 5 for the. constitutive 
model presented in Appendix 1 The physical vessel dimensions were 7-mm outer 
diameter and 4-mm inner diameter as to approximate the physiological case of a 
femoral artery. Figures 24a and. 24b respectively present the lateral and axial 
displacement fields. Maximum motion opcured at the lumen- Figures 24c to 24f 
present the Aij components of Eq. 7;. which are. respectively the lateral strain, the 
lateral shear; the axial shear and the axial strain. 



The elastogram,. i.e. the. map of the axial strain distribution, is generally 
presented as a grayscaie.; image. Because in conventional elastography, external 

5 . 

y-axis compression is applied and also because 1D y-axis motion is assessed, Ayy 
is expected £ 0. Traditionally, smaller strain amplitude values are associated with 



5 Notice that 6 % intraluminal dilation ts equivalent to 3 % compression of the intraluminal wail. 
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harder , regions and are printed. In bfack; equivalents, . higher strain amplitude . 
values, are associated with softer regions and are printed in white. However in . 
NIVE, as it can be observed in Fig. 24f, dilation can also be detected (Ayy £ 0) in 
the elastogram. in an elastographic sense, the dilation, regions can be 
misinterpreted as soft tissue. Indeed; in Fig. 24f, two harder zones (Ayy < 0} likely * 
seem to be identified at 12 and 6 o'clock. Because, for the conditions simulated, 
the vessel wall is homogeneous, such a phenomenon is . referred to as hardening 
artefact. Inversely, two softer zones (Ayy £ 0) at 3 and 9 o'clock seem also to be 
identified; they are referred to as soften tag artefacts. Such motion artefacts stem 
from the fact that motion occurs radially and is observed in Cartesian coordinates; 
this can constitute a potential limitation to lion-invasive vascular elastography. 

To go further with, this observation, let's consider the motion parameters in 
their natural polar coordinate system, or material coordinate system. In Fig. 25a, is 
presented the radial displacement field , computed from the lateral and axial 
displacement fields (Figs 24a and 24b, respectively); This radial displacement field 
is also presented in a polar (r,<p) coordinate system (Fig. 25b). The gradient of the 
latter displacement field thus provides. the radial strain (Fig. 25c). In Fig. 25d is 
shown a plot of the radial strain at <p = w. One can observe the mdnotonic profile of 
this plot, bein£ maximal at the lumen and minimal at the outer side of the vessel; 
such a phenomenon is a consequence qf the boundary conditions and is known, in 
the literature, as the strain decay. Finally in Fig. 25e, the radial strain is reported 
back in the (x,y) coordinate system. Regardless of the strain decay phenomenon, 
one can observe that no specific hard or soft region is identified. Fig. 25e thus 
illustrates a strain profile that is quite representative of a. homogenous vessel wall 
behavior, 

A new parameter for tissue characterization 
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Hopefully, elastograms such as the one in Fig. 25e are required to : 
appropriately characterize .the vessel wall. However, to prevent signal ■ 
decorrelation in NiVE, motiqn is to be studied in the transducer coordinate system; 
that is the (x,y)-Cartesian coordinates. Accordingly, elastograms are expected to : 
be as artefabtual as the one In Fig. 241 However, taking a close look at all the ] 
motion parameters in Fig. 24,. one can observe that A** and Ayy are 
complementary, while Axy = Ayx. In this paper, the Von Mises (VM) coefficient is I 
proposed as a parameter to, characterize the vessel wall. Indeed, VM is ; 
independent of the coordinate; system arid is a combination of the four j 
10 displacement gradient components (Ay) that is mathematically expressed as: 

(9) 

In Fig. 26, a comparison between 4 and the radial strain (Fig. 15e) is ; 

presented. Qualitatively, both parameters are equivalent. Fig. 26c shows the plots i 

for the radial strain ( — ) and E; (--) at x ~ 0. Analyzing those curves, one can ; 

observe , that the VM coefficient likely improves the contrast between higher and : 

lower strains: while the profile remains the Same. Moreover, regardless of the strain . 

20 decay, £ (as well as the radial strain) is interestingly free of hardening or softening : 

artefacts. It i$ thus allowed believing that such a. parameter could be very suitable i 

to nonrinvasively characterize the vessel wall. To corroborate such an assumption, ■ 

a more complex geometry, that is $ heterogeneous vessel wall, is considered in ] 
the next section. 

Characterization of an heterogeneous vessel wall 
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The case of a pathologic yesse! wall ?s now investigated, that is when the . 
lumen is narrowed by the presence of a hard plaque at the intima. It is assumed : 
that the plaque is 10 times stjffer than the normal vascular tissue, a lipid pool 10 \ 
times softer than the normal vascular tissue immediately surrounds it. The 
geometry is schematically presented in Fijj. 27a. to have approximately the same : 
range of strains as in the homogenous case study, the intraluminal pressure ; 
gradient was . set at 40 rhmHg. (intraluminal dilation close to . 6 %). The. outer 
diameter of the vessel was again set at 7 mm. In Figs 27b to 27e, the- four Ag ; 
tensor components are presented. One can appreciate the very complex patterns 
10 of those parameters. The elastogram (i.e; t the axial strain map), presented in Fig. 
27e (equivalent^ the lateral strain map of Fig. 27b), shows a region of relatively 
high strain values; this corresponds to the location of the lipid pool as it cm be 
expected; Nevertheless, the elastpgram does not undoubtedly allow differentiating 
between the plaque, the lipid pool and the normal vascular tissue. 

In Fig. 28b, the elastogram computed using the VM parameter {§) is 
presented. Obviously, three. main regions are identified: the hard plaque with lower 
strain values, :the lipid pool; with higher strain values and the normal tissue with . 
moderate strain values. This elastogram can qualitatively be compared with the 
20 radial strain map presented in Fig, 28a. For a quantitative comparison, Fig. 18c 
shows plots of the radial strain and § at x = 0. In both of them, the three specific 
regions can be clearly differentiated.. As it was also observed in the homogeneous 
case study, it is interesting , to. point out that \ substantially improves the contrast 
between soft and hard tissues. : 

These observations are interesting and allow believing in the potential 
application of NIVE. Nevertheless, some artefacts still remain, specifically at the 
intersections: between haFd, soft and. normal tissues. This is, unfortunately, quite 
expectable since those are sites of maximal stress concentrations. Such a 
30 problem couljd have been circumvented if the Young'&modulus, instead of a strain 
parameter, was computed in ultrasound eiastography. Another problem, which 
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could, arise in the real-life, situation, concerns the limited . lateral resolution of 
current ultrasound systems. Because of that, variances in the lateral motion . 
estimates and Axy) may in practice be large. This may constitute a potential 
limitation to the assessment of In the -next sections, . we present results from 
simulated RF data* All elasfograrris were compMted with the Lagrangian speckle 
model estimator (LSME). 

RESULTS 

10 

APPLICATION OF THE LSME TO SIMULATED RF ECHO SIGNALS 

The dynamic imagei-fcfmation model of Eq< 5 was used to simulate 
sequences of R F images. As illustrated in. Fig. 21, the vessel section (the 
scattering function Z(x,yj) was modeled as a 2D white Gaussian noise (number of 
scatters in a resolution cell. » 5); The lumen and the tissue surrounding the vessel 
were assumed tp.be respectively 2-5 anil 1.67 times less echoic than the wall. 
Z(x,y) was low-pass filtered before considering motion; the post-motion scattering 
function was expressed as Z tT ., (x, y) . The radial displacement field used to mimic 

20 motion was? obtained from the analytical bfomechanical model presented in 
Appendix I (normal cylindrical vessel), and Ansys simulations (obstructed vessel 
with a plaque). It is important to remember that such a motion is normally Induced 
by the cardiac pulsation. The imaging system PSF (H(x,y)) was defined as a 2D 
Gaussian wavelet, modeling a 10 MHz transducer with a 60% bandwidth at ~3dB 
and a beam width. (width at half maximum) of 0.6 mm. The simulated RF images 
(pre- and post-tissue motion images) vvere used as inputs by the LSME to 
compute the. four displacement gradient. components. Are presented, in the next 
two sections, results for the homogenous snd the. heterogeneous vessel wall 
studies, respectively. 
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Motion estimation for a homogenous vessel wall 

As defined in section 2.4.1, . a vessel of 7-rnm outer diameter and 4-mrn 
inner diameter was simulated. To take into account the surrounding tissue, the 
echo field was 8 mm x 8 mm; that is 3000 samples axially * 150 RF lines laterally. 
To apply the LSME, a measurement-window of 533 p.m x 1 066 j*m (200 samples 
x 20 RF lines) was chosen. As illustrated in Fig. 20, the motion was estimated for 
windows W mn with 88 % and 90 % axial and lateral overlap, respectively, 

10 The windowed data of the pre and post-tissue-motion images were pre- 

processed to compensate for tissue translation using 2D correlation. This left the 
LT matrix and eventually a small residual translation motion (61 and 64) to 
estimate. In this, paper, the ..translation motion parameters 9t and '64 were not 
considered. The four displacement gradient components, as computed with the 
LSME, are presented in Figs 29a to 29d. 

Figs 29a and 29b . present the maps of the lateral strain and shear, 
respectively. As. it can be expected, those motion parameters show , large 
variances. That is due to the limited lateral image resolution (53 pm per RF line, 
20 for this study) used to simulate the performance of current ultrasound systems. 
Nevertheless,; the lateral strain values are observed in a range, on average, that 
includes the theoretical results (Fig 14c): Figs 29c and 29d present; maps of the 
axial shear and strain, respectively. The elastogram (Fig. 29d) is quantitatively 
similar to the theoretical one (Fig. 24f), showing strains in the interval of [-3 %, 3 
%]. The pattern of the axial -shear is also quantitatively similar to the theoretical 
one (Fig, 24e). 

Are presented in Fig. 30, comparisons between the simulated and 
theoretical VM parameters For a quantitative comparison, Fig. 30b presents an 
30 average of three curves around x = 0, and at Fig. 30d the. curve corresponds to x = 
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0. Notice the presence of strain decay in each curve. It is also important to 
emphasize that close to the lumen, motion estimates, may be, in some instances, ; 
inaccurate. That is explained by the fact that in such locations, the measurement- 
windows may overlap the vessel wall (where pre- and post-tissue-mbtion signals 
are expected to be coherent) and the lumen (where pre- and post-tissue-motion 
signals are expected to be. uncorreiated). As a consequence of that such motion 
estimates may not be reliable; Because one is now aware of some of the eventual j 
difficulties that could be encountered and,, in the. light of the promising results of 
Fig, 30, it allows believing .in the potential of N1VE. In the next section, the 
10 heterogeneous vessel wail is investigated. 

Motion estimation for a hetemgeneous vessel '.wall 

The vessel wall investigated is the same as in section 2 A3. The data ; 
acquisition system (i.e. the PSF). and. the parameters for motion assessment (W mfh ; 
overlap, etc.) are the same as in the section 3.1.1 above. In Fig. 31, the results for ) 
the four displacement gradient components are presented. Figs 31a and 31b 
present the -maps for the lateral strain and shear, respectively! As for the 
homogenous: case, the variances for those motion estimates are large; they do not ; 

20 clearly allow differentiating between hard and soft structures. Nevertheless, a 
reasonable range of strain values is, on average, observed when comparisons are : 
made with . Fig. 27. The simulated elastogram (Fig, 31 d) shows a slight 
underestimation of the axial strain valued when compared to the theoretical one 
(Fig. 27e); strain values in the interval of [-2,5 %, 2,5 %] are observed instead of [~ 
6 %, 6 %]. That is, at feast partially, due to the fact that high strain areas at the ! 
intersections; between hard,, soft and normal tissues are more prominent in the - 
theoretical ^lastogram thari in the simulated one. The axial shear is also = 
qualitatively and quantitatively similar to the theoretical one (Fig.. 27d). In Fig. 32a, 
is shown the.; composite elastograrp (£). Even with the presence of large variances 

30 on the estimates of the lateral motion parameters, § allows identifying the hard ' 
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plaque and thfe soft lipid pool* For a quantitative validation, Fig. 32b shows an 
average of three curves around x =' 0. Clearly the hard plaque (strain values less \ 
than 1 %) and the lipid pool (strain values around 4 %) are detected. Because of ': 
the strain decay phenomenon, the normal tissue presents strain values standing 
from 1 % up to 2 %. For comparison purpose, the same results are shown for the 
theoretical \ in Figs 32c and $2d. . 

DISCUSSION 

10 ; • : 

The vascular tissue is made of elastin, collagen and smooth muscle cells; its 
mechanical properties are thus very: complex. Pathological conditions of such a 
tissue often: induce changes in the vessel, wall elasticity. For example, plaque ■ 
deposit stiffens the vascular jwall and theni counteracts its dilation under systolic 
blood pressure. On the other hand, peripheral aneurysmal rupture of the aorta has : 
been linked to changes in the proportion and integrity of collagen and elastin* 
Because of that, vessel diameters as well as their variations during the cardiac 
cycle were considered as indices to characterize vessel pathologies. Hence, from • 
a* biomechanfcal perspective,, it can be understood that to better appreciate the . 

20 dynamics of the arterial wall and its pathologies, a more detailed description of the 
mechanical and elastic properties of the arteries is required. 

Elasticity is a very suitable parameter to describe vessel wall function. It may • 
be quantified; with compliance, a measurement that expresses the relative change : 
in vessel crgss-sectional area as|a function of time, However, a rrtare complete ; 
method to outline the elastic properties of the vessel wall seems to be the 
elastography; indeed, sufch a technology aims to provide images of the elastic 
properties of biological soft tissues, 

30 This, paper investigated the feasibility of NIVE. It was shown that a major 

difficulty with this approach ;stems from the fact that motion oqcurs radially within • 
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the. vessel wall while the ultrasound beam propagates axialfy (linear array 
transducers were considered in the present study). Regarding that, a 2D tissue- 
motion estimator was required. Moreover,: in such a situation, the motion 
parameters are expected to be very difficult to interpret One major challenge with 
NIVE then concerns the interpretation of the motion estimates to appropriately 
characterize the vascular tissue, 

ABOUT THE FORWARD PROBLEM 

As a first step, the forward problem (FP) in NIVE was; addressed. 
Biomechanical simulations of the vascular wall properties were performed and 
motion parameters were investigated. The FP allowed to better understanding 
vascular tissue dynamics in NIVE. Indeed, the analytical solution for a 
homogeneous vascular section subjectedrto an intraluminal pressure was derived. • 
Ansys simulations were also performed to;. address the more complex, situiation of a 
heterogeneous plaque: Because tissue motion occurs radially, within the vessel I 
wall while the motion is assgssed : in.(x,y)-Cartesia.n coordinates, the elastogram • 
(Ayy) was shown to be subjected to: hardening and softening artefacts; 

To overcome the problems related to those mechanical artefacts, a new 
characterization parameter (£) was proposed. Such a parameter, also known in 
the literature as the Von Mises (VM) coefficient, uses the full 2D-strain tensor to ; 
provide a measure free from hardening! and softening artefacts. % likely allows \ 
characterizing the vessel waif in its natural coordinate system. Indeed, such an 
elastogram presents some similarities with the radial strain map, and even shows 
up the very, well known strain decay as it should be expected. For the case of the I 
simulated pathological vessel, the VM parameter successfully allowed \ 
differentiating between hard, .soft and normal vascular tissues, NIVE should then 
be considered, hi some instances, as. a promising tool for non-invasiveiy; 
investigating vessel walls. 
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ABOUT THE APPLICATION OF THE LSME 

A 2D-tiss:ue motion estimator, known in the literature as the Lagrangian . 
speckle estimator (LSME), Was used to. assess motion. The LSME is; a 2D niodel- ; 
based estimator that allows computing the full 2D-strain tensor. The approach ; 
relies oh a .tissue-motion mpdel and on a dynamic image-formation model. It : 
assumes that the image formation can be modeled as a linear space-invariant 
operation on a scattering, function,. Indeed, the simulated RF images result from 
the 2D-convo;Iution operation between th£ point-spread function of an ultrasound : 
10 system and the acoustic impedance functions, which were modeled as a white 
Gaussian noise. 

The same data sets as far the FP were investigated- As well, for the : : 
homogeneous and heterogeneous, vessels, the results seemed promising. The = 
Von Mises characterization parameter (ty showed great potential to differentiate 
hard , and soft tissues from the normal vascular one. However, the main problem 
with the LSME remains the assessment of the lateral motion parameters. As a \ 
matter of fact, due to the relatively, poor: lateral resolution of RF images, the j 
variance for those motion estimates is large, it is also important to emphasize that 
20 the 2D motion estimator, presented in this paper, assumes that motion occurs in \ 
plane strain condition. However ^practice, because of blood pressure excitation, • 
the vessel wall is also expected to be under distension stress. As a consequence ■ 
to that, out-;of-plane motion will occur in ;the dynamics between the,pre- and the '" 
post-motion RF images; this is another factor that may increase the variance of ; 
motion estimates. To at least partially counteract out-of-plane motion, one could ; 
ensure that a small intraluminal pressure: gradient is induced, This could be done . 
usiing an ECG gating along with the ultrasound acquisition system. The 
a Lagrangian":filter was also shown.to.be effective to compensate for such a source \ 
of decorrelation noise. 
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Another potential issue to vascular elastography. is that vessel walls are 
subjected to residual stresses. Indeed, it was demonstrated that the magnitudes of 
such residual ^stresses (equlvalentjy residual strains) may be relatively Significant. 
Accordingly, vascular elastograms are expected to be fully valid whenever residual 
strains are much smaller than strains as induced by the intraluminal pressure 
gradient. Provided that, in NIVE, ejastograms compute maps of relative stiffness, 
arid provided ; that tissue motion is Assessed, over small regions* of interest where 
linearity is assumed, . residual stresses may not have significant effects. 
Furthermore, another potential solution may consist in setting gri optimal pressure 
10 gradient such as the tissue strains are much larger than potential residual strains, 
but are small enough to be in/the limits of Shear elasticity. 

CONCLUSION 

An approach to nornnvasively characterize vessel walls, labeled as NIVE, 
has been introduced. Whereas endova^cular elastography is invasive and its 
clinical application restricted to a complementary tool to assist 1VUS echograms in 
pre-operativfc: lesion assessments and to plan endovascular therapy, NIVE method 
20 may be of value to investigate vascular tissue properties of patients under 
medication, and for post-surgioa! follow-iiip v Additional to that, NIVE could be of 
interest to investigate small vessels where catheterization and IVUS echography is 
not possible. While the use of a balloon catheter is not required in NIVE, such an 
option can be useful if a more precise. dilation increment is needed. In this case, 
however, the technique would be invasive. Regarding the forward problem, the 
Von Mises coefficient has been proposed as a new characterization parameter. 
Numerical simulations showed the. potential of the VM parameter to differentiate 
between hard, soft and normal vascular tissues. The lateral resolution of current 
ultrasound systems is, however, a limitation to such an approach. 
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FIGURE CAPTIONS 

Figure 20.- 'Simplified illustration of NIVE (non-invasive vascular elastography) 
acquisition system, showing that the tissue motion occurs radiMy, whereas the 
linear-array transducer scans, in the Cartesian coordinate system.; Processing 
motion estimation requires 2D segmentation of the images; the Wm* represent the 
data windows that are used by the Lagrangisin $peckle model estimator. 

Figure 21 Schematic implementation of the dynamic image-formation model. For 
10 this example, thfc displacement field was qomputed from the analytical solution for 
a pressurized thick-wall cylindrical bipod vessel embedded in an elastic coaxial 
cylindrical medium (Appendix I); The parameters on this figure refer, to Eq. 5. The 
dynamic image-formation model was . also applied to the case of the 
heterogeneous vessel. In this cas£, the displacement field was determined with 
Ansys simulations. 

Figure 22.- Block-diagram showing, the implementation of the Lagrangian speckle 
motion model estimator (Eq. 6). The algorithm searches the tissUe motion vector \|r 
that best matches l 0 and is the initial guess that is required to start the 
20 iterative process. In this study, a 2 % axial compression was : assumed for the 
initial guess: 

Figure 23x illustration of the difference between the observer's (x,y) and the 
material (r,cp) coordinate systems. The flatter, also known as the Lagrangian 
coordinate system, considers each portion of a continuum as if it w£re a particle, 
for which the :trajectory is described as a f unctio n of time. 

Figure 24.-;Motion parameters for a pressurized thick-wall cylindrical; blood vessel, 
embedded in an elastic infinite medium. Rigs 24a and 24b present respectively the 
30 lateral and axial displacement fields; the colorba.rs express the displacements in 
H.m. Figs 24c to 24f show the Aj| components of Eq. 3, which are respectively the 
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lateral strain (Axx), the lateral shear (A xy ), the axial shear (Ayx) and the axial strain 
(Ayy);. the qolarbars express the strain in percent 

Figure 25.- Radial strain and strain decay for. a homogeneous vessel wall a) 
Radial displacement computed frprh the lateral and axial displacement fields; b) 
radial displacement field in a polar (r^) coordinate system; c) radial strain 
computed from the gradient of 25b; d) radial strain at <p » n showing the strain 
decay; e) radial strain reported hack in the (x,y)-coordinate system. 

10 Figure 26.- Comparison between the radial strain and the Von Mises parameter 
(£) for a homogeneous vessel, wall a) Radial strain in the (x,y)-coordinate system; 
b) : Map of £; c) curves from the radial strain (— ) and £ (— ) at x « 0. 

Figure 27.- Motion parameters for an heterogeneous vessel wall; a) Simplified 
representation of the geometry; b) to. e) show the A9. components of Eq. 3, which 
are respectively the lateral strain (Axx), the lateral shear (Axy), the. axial shear (Ayx), 
and the axial; strain (A yy ); they were computed from the gradients: of the lateral and 
axial displacement fields. The;colorbars express the strain in percent. : 

20 Figure 28*- Comparison. between the radial strain and the Von .Mises parameter 
a) Map of the <x f y)~radM strain;, b) map of §; c) curves of a) (— ) and b) (— ) at 
x = (X 

Figure 29.- Motion parameters as .computed with the Lagrangian speckle model 
estimator for a homogenous vessel wall; a) to d) Ay components of Eq. 3, which 
are respectively the lateral/strain (A xx ) f the : lateral shear (A xy ), the axial shear (Ayx), 
and the axial; ^train (Ayy). 

Figure 10,- Comparison between simulated and theoretical Von Mises parameters . 
30 (§) for a homogeneous vessel . w0. a) Map of . the simulated £; b) curve of the 
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simulated £ at.x = 0; c) map of the theoretical §; d) curve of the theoretical § at x = 
0. 

Figure 11.- Motion parameters as /estimated with the Lagrangian speckle model 
estimator foran heterogeneous vessel wall; a), to d) Ag components of Eq. 3, which 
are respectively the lateral strain (Axx) t . thefateral shear (A X y), the axial shear (Ayx), 
and the axial strain (Ayy). 

Figure 12.- Comparison between simulated and theoretical Von Mfees parameters 
(4) for an heterogeneous vessel w^lf. a) Map of the simulated £;; b) curve of the 
simulated £ at x = Q; c) map of the theoretical ^; d) curve of the theoretical £ at x - 
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APPENDIX I 

Thick-wall cylindrical blood vessel embedded in a. finite (or infinite) medium 

An exact solution of a pressurized thick-wall cylindrical blood vessel of inner 
and outer radii Rj and R 0 . respectively, embedded in an elastic coaxial cylindrical 
medium of radius Re f can be found in linear elasticity. For the mathematical 
formulation, the cylindrical and Cartesian unit base vectors and their associated 
physical coordinates are .noted (e r ,e^ t e]) $ (e x ,e yt e s ) and (r, <p, 2), (x, y, z), 

10 respectively. In this mechanical problem, the Cartesian vessel wall displacement 
field and strain components are of interest. The assumption of plane strain was 
made (in the (r t <p)-plane) because the vessel length is at least of thd same, order 
of magnitude as its radial dimension. 

The twp elastic media are. assumed to be incompressible arid isotropic and 
are described by the. constitutive laws; 

[o-Y*-^ (A.1) 

20 where the superscript 'm' denotes ;the considered medium (m=1 for the vascular 
wall for which &£r£'jt 0 -and 0<q>^27c l -and m=2 for the embedded tissue for 

which R„<r£Re and Q£<p£2n); The parameters [cr] my and {ef m) are the stress 
and strain tensors, [/] is the identity matrix, E im) are the Young's moduli, and 
p im) are the tagrangian multipliers resisting from the incompressibility of the 
materials [38] given by the. following kinematics constraints: 

V..w ( ! w) -0 ' ' . (A.2) 
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where u (m) are unknown. displacement vectors! If gravity and inertia! forces are 
neglected, the conditions of local equilibrium are V. [&f m) ==0, or in terms of the ; 
displacement vectors: 



; 3 



The displacement fields u in \ as weli> as the stress tensors [af m) 
the following boundary conditions. The blood pressure P b is uniform, 



10 



at r = R| 



(A.3) 



must satisfy 



(A.4) 



no stresses are applied on the external surface of the surrounded, tissue, 

[df^e^O at r = R e (A,5) 



arid at the interface between the vessel wall and the surrounded tissue, we must 
have the continuity of displacement fields and equality of the stress vectors, 



20 



at 

af ; 



r = R 0 
r= R 0 



(A.6) 
(A.7) 



Due to the symmetry, of the mechanical problem and. because the two 
medial are incompressible^ the. displacement solutions are: 



r 



(A.8) 
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where K is a constant. The equilibrium conditions (Eq. A.3) are satisfied only if the 
two Lagrangian terms jc? (,) and ,p i2) are constant. At the end, the three unknowns of 
the problem (K. p il K p m ) are found by using the boundary conditions, given by 
Eqs. (A,4), (A.5) and (A;7). Hence, for the value of K, we have; 



2 *[ [R? R]); ■ \# Rl) 



(A.9) 



In the 'particular case of a thick-wall cylindrical blood vessel embedded in an 
infinite medium, the new value of K is obtained from the limit = Iim K , and is 

10 given by: .• ■ < 



(A. 10) 



So, the needed Cartesian components of the displacement vectors and strain 
tensors are: . 



(A.11) 
(A-12) 



20 



2 3y dx. 



M (^:+:/) 2 



(A.13) 
(A.14) 
(A.15) 
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FOOTNOTES 

1 This is the Matlab notation fpr augmented vector (;) and matrix vectorisatioh .(:). 
Hence, \j/ is a 6 x 1 vector built from the 2 x 1 Tr vector and the 4x1 vectorisation 
of LT. 



2 Notice that 6 % intraluminal dilation is equivalent to 3 % compression of the 
intraluminal wall. 
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CLAiMS 

1. A method for characterising the mechanical properties of tissues in smalt 
vessels using the Non-Invasive Ultrasound Elastography (N1VE), the method 
comprising the steps of : 

acquiring pre- and post- deformation images; 

using the Lagrangian speckle tissue model estimator to assess the 2~D 
stain tensor; 

computing the Von Mises elastrograms; and 
10 modelising the tissue mechanical characteristics. 

2. An apparatus for investigating the feasability of N1VE for characterising tissues 
in small vessels. 

3. A model-based method devoted to endovascular elastography, the method 
comprising the steps of : 

performing a modeling in performing the sub-steps of : 

A1) modeling the polar static-image formation; 

A2) modeling the tissue-motion; 
2 0 A3) modeling the polar dynamic-image-formation; 

estimating tissue motion in performing the sub-steps of : 

B1) assuming speckle is a material property; 

B2) providing LSI (Lagrangian speckle images); 

B3) using the Lagrangian Speckle Mode! Estimator (LSME) to assess 
tissue motion; 

B4) solving a non linear minimization problem; 

B5) using the optical flow equations; 
thereby providing an accurate tissue motion estimate for providing an 
elastogram of the tissue. 

30 
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4. An apparatus for characterising the mechanical properties of tissues, the 
apparatus comprising : 

a CVIS ultrasound scanner; 

a 30-4Q MHz mechanical rotating single element transducer; 
an oscilloscope; 

a pressuring system; and 
a computer. 

5. A method for detecting the hard plaque in the blood vessels, the method 
comprising the steps of : 

collecting a sequence of RF images while incrementaiy adjusting the 
intraluminal static pressure;and 

characterising the the hard plaque from the normal vascular tissue. 

6. An apparatus for detecting the hard plaque in the blood vessels, the apparatus 
comprising : 

a CVIS ultrasound scanner; 

a 30 MHz mechanical rotating single element transducer; 

an oscilloscope; 

a pressuring system; and 

a computer. 
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Fig. 2) Schematic illustration of the vascular flow 

phantom. 



Cft. 02457171 2004-02-09 




Fig. 3) Schematic representation of the moulds that were 
used to manufacture the double-layer PVA-C vessel-mimicking 
phantom. The simulated vessel had a 1 . 5 -mm lumen diameter, 
1 . 5 -mm wall thickness (approximately 0.75-mm for each 
layer), and 5 2 -mm length. 



Transducer 




Fig. 4) Simplified illustration of NIVE acquisition system, 
showing that the tissue motion occurs radially, whereas the 
transducer scans in the Cartesian coordinate system. The 
approach assumes a constant deformation for each W^. 



Pre-motion image 



■ i 



Post-motion image 



Map of the correlation coefficient 



Von Mises parameter 




Fig. 5) a) and b) show pre- and post-motion B-mode images 
(RF images were processed to achieve tissue motion) ; c) 

presents a map of the correlation coefficient between the 

" pre- and post-motion RF images; d) shows the composite 
elastogram (*) as computed by the LSME (the color bar 
indicates the strain in %); e) presents a zoomed axial 

portion of d, showing the inner layer (softer) with higher 
strain values, and the outer layer (harder) with lower 
strain values; f) plot of an average of 20 axial lines 
chosen in the middle of %. 
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FIG. 6 




FIG. 7 




FIG. 8 



Cross-sectional vessel wall reconstruction 
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FIG. 15 
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Fig. 16) Experimental set-up for intravascular 
elastographic investigations. 
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a) 



the plaque from the healthy vascular tissue. 
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Fiq 18)- a) to j) are elastograms computed for consecutive 
Increasing physiologic fluid pressure levels Lateral and 
Sia! vilue? are dimensions, while the colorbars gxve the 

strain in percent. 




Fig. 19)- a) Histological section, showing a very small 
atherosclerotic plaque at 3 o'clock; b) Elastogram, showing 
very significant qualitative and quantitative similarities 

with histology. 
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Figure 21 
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Heterogeneous vessel wall geometry 
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Figure 27 
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